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PREFACE 


The initial development of this analysis was conducted under Corporate 
sponsored independent research and development funding at United Technologies 
Research Center. Extensive refinements to the analysis were made under spon- 
sorship of the Langley Research Center of the National Aeronautics and Space 
Administration and the U. S. Army Air Mobility Research and Development Labo- 
ratory, Langley Directorate as part of Contract NASI-IO 96 O. The purpose of 
this contract was to perform an investigation of a bearingless helicopter 
rotor concept having a composite primary structure. A coopanlon report to 
the one herein is NASA CR- 2637 which presents the results of graphite/epoxy 
fatigue tests, wind tunnel experiments, correlation studies, and a preliminary 
design of a full scale helicopter rotor, also obtained under this contract. 

The correlation studies were used to provide limited validation of the herein 
described analysis, and the preliminary full scale design study involved sub- 
stantial use of the analysis to investigate the structural dynamics and aero- 
elastics of the design. 

The author wishes to express his appreciation to Mr. C. E. Swindlehurst 
and Dr. W. F. White for their encouragement and useful comments. 
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AEROELASTIC ANALYSIS FOR HELICOPTER ROTOR BLADES 
WITH TIME-VARIABLE, NONLINEAR STRUCTURAL TWIST 
AND MULTIPLE STRUCTURAL REDUNDANCY - MATHEMATICAL 
DERIVATION AND PROGRAM USER’S MANUAL 


Richard L. Bielawa 
United Technologies Research Center 


SUMMARY 


The differential equations of motion for the lateral and torsional 
deformations of a nonlinearly twisted rotor blade in steady fli^t conditions 
together with those additional aeroelastic featvires germane to composite 
bearingless rotors are derived. The differential equations are formulated in 
teims of uncoupled (zero pitch and twist) vibratory modes with exact coupling 
effects due to finite, time variable blade pitch and, to second order, twist. 
Also presented are derivations of the fully coupled inertia and aerodynamic 
load distributions, automatic pitch change coupling effects, structural redun- 
dancy characteristics of the composite bearingless rotor flexbeam - torque tube 
system in bending and torsion, and a description of the linearized equations 
appropriate for eigensolution analyses. Three appendices are included pre- 
senting material appropriate to the digital computer program implementation 
of the analysis, program g400. 



INTRODUCTION 


The composite bearingless rotor employs for its primary structural 
element a spar fashioned from radially aligned uniaxial high strength fibers 
(carbon, boron, etc.) in an epoxy matrix. The transverse shear modulus of 
such a spar is sufficiently low to produce a torsionally flexible member which 
when installed over a finite length (designated the flexbeam ) replaces the 
feathering bearings normally used for blade pitch control. Figure 1 shows the 
schematic of a typical rotor employing such a spar. Blade pitch control is 
achieved by elastically twisting the inboard portion of the spar; the moment 
applied to the blade from the push-rod is transmitted through the aerodynamic 
shell or torque tube which, like the outer portion of the blade, is relatively 
stiff in torsion. 

These structural features produce unique aeroelastic characteristics not 
readily amenable to conventional analysis; First, the primary structural 
member, the spar and included flexbeam is subjected to highly nonlinear and 
time varying structural twist which, over the flexbeam span, can approach 
(equivalent) total span twist angles of +90“ . Second, although the torque 
tube and flexbeam have contrasting specialized structural functions (torsion 
and bending load transmissibilities , respectively) each nonetheless exhibits 
significant amounts of both types of load transmissibility . Hence, the torque 
tube - flexbeam system comprises a doubly redundant structure: redundancy both 

in torsion and in bending. Third, the simple torque tube shown in Fig. 1 is 
that of a "cantilevered" configuratiai wherein the inboard end is supported 
both in shear and torque solely by the push-rod. Such a configuration pro- . 
duces not only a soft blade torsion system, but a high degree of pitch-flap 
coupling . 

Because of these resulting unique aeroelastic characteristics, none of 
the various "comprehensive" (nonlinear) rotor aeroelastic analyses currently 
available (e.g.. Refs. 1, 2, and 3) can be Justifiably applied to the bearing- 
less rotor. Most of the difficulties encountered with these analyses are due, 
in one form or another, to inappropriate assumptions made with regard to struc- 
tural twist. Compared on this basis, the various analyses generally fall into 
either of two categories. In the first category, the structural twist is 
assumed to be linear, small, and temporally constant; such an assumption is 
clearly at variance with the aforementioned twist characteristics of bearing- 
less rotors. In the second category, large nonlinear structural twists are 
incorporated by the use of "coupled modes" wherein the blade bending elasticity 
is defined by the use of the normal (vibrational) modes of an arbitrarily 
pretwisted blade at some nominal fixed collective angle. Unfortunately, such 
an approach is rigorous only in hover wherein the control angle is constant and 
equal to the nominal collective angle. Generally, the use of coupled modal 
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analyses for typical forward flight cases involving substantial cyclic angles 
leads to an obvious contradiction of the inherent assumptions upon which 
coupled modes are predicted. To rectify this conceptual deficiency, coupled 
mode analyses must incorporate, in some fashion, the mode shape and modal 
frequency variations over the appropriate pitch angle range as defined by the 
given flight case. This invariably leads either to extreme complexity and 
resulting computer storage and run-time requirements, or to various simplifying 
approximations to the variations of the mode shapes and frequencies with pitch 
angle. Either of these consequences, however, tend to nullify the advantages 
claimed for "coupled mode" analyses. In Ref. 3 is presented a consistent sys- 
tematic development of such a coupled mode analyses which rigorously addresses 
itself to the problem arising from variable pitch angle. However, even for 
this exemplary analysis, the time variable (perturbational) pitch angle is 
assumed to be small, constant with span and, hence, inappropriate to the 
requirements of con^osite bearingless rotors. Because the use of coupled 
modes offered no clear advantage and because a valid, relatively simple yet 
practical means of incorporating nonlinear structural twist was found, a 
coupled modes approach was discarded in favor of one based upon more conven- 
tional uncoupled modes. Finally, in all cases the available analyses were 
found to be totally incapable of providing an analysis of the structural 
redundancy of the torque tube - flexbeam system or of a cantilever configured 
torque tube. 

The aeroelastic analysis described herein is a multi-purpose computer 
program characterized by rigorous modelings of nonlinear and time varying 
structural twist, and of the redundant load carrying features of the bearing- 
less rotor. Although developed in response to the requirements of composite 
bearingless rotors , the dynamic equations are sufficiently general for valid 
application to all conventional rotor systems: articulated, semiart iculated, 

teetering and hingeless. The computer program implementing these equations 
presently assumes a fixed hub, and hence, obtains solutions of the dynamic 
equation for only one blade. It cannot, therefore, be used for teetering 
rotor systems which inherently require a two-bladed implementation of these 
equations. The differential equations of blade beam bending (flatwise and 
edgewise) and torsion are solved using a Galerkin procedure wherein the normal 
"uncoupled mode" shapes and spanwise derivatives of the blade pitch angle and 
the nonlinear twist are appropriately coinbined to describe the "coupled" blade 
deflections. The general approach to the development of the aeroelastis anal- 
ysis closely parallels and draws upon that used and reported in Ref. 1. Two 
types of solution are available: eigensolutions of various linearized equa- 

tion sets for coupled frequency and/or stability analysis purposes, and time- 
history solutions of the complete nonlinear equations for harmonic analysis 
and/or transient aeroelastic response calculation purposes. The aerodynamic 
description includes the use of predetermined static airfoil data, constant or 
variable inflow (vorticity induced and/or momentum derived) and unsteady dynam- 
ic stall data. 


This document presents the more salient details of this aeroelastlc 
analysis. Specific Items to be described are: 1) the principal Inherent 

assumptions, 2) the coordinate transformation used to Introduce pitch and 
twist effects, 3) application of the Galerkln procedure to obtain the basic 
modal equations, 4) descriptions of the required load distributions, Intro- 
duction of pitch -flat /edge coupling, 6) simulation of the '^wobble mode," 

7) Inclusion of the redundant load carrying features of the flexbeam - torque 
tube assembly of the bearingless rotor, and 8) the detailed statement of the 
linearized equations used in the elgensolutlon. 

Ilhese elg^t sections essentially present the derivation of equations of 
motion whose solutions can then be Implemented using a variety of numerical 
techniques. In Appendix I Is presented a brief description of the quadrature 
techniques employed In the existing digital computer program Implementation of 
these equations, program g 400. Appendix II contains a detailed description of 
the Input required to run this digital computer progrcun, and Appendix III pre- 
sents a comprehensive description of the cooqputer generated output. 
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LIST OF SYMBOLS 


SgO 

A 
A 


B 


c 


Cd 


Acd^ 


H 


C 

C, D, 
^ld 

Cj 


D 

e 


®A 


free stream sonic velocity, 
area of blade section, 
inertia matrix. 

general submatrices of partitioned inertia matrix without and 
with the effects of coupling due to automatic pitch change, 
respectively. 

tip loss factor, used in momentum inflow equations and in 
approximate simulation of three-dimensional tip effects. 

blade or section chord, as appropriate. 

section drag coefficient 

increment to sectibn drag coefficient to account for surface 
roughness. 

section lift coefficient. 

section pitching moment coefficient about quarter chord. 

compression at an arbitrary spanwise station. 

integration constants from solution for flexbeam torsional 
element differential equation. 

rolling and pitching moment coefficients, respectively, of the 
rotor; mcments nondimensionalized by (p'lrE^Cp') • 

coefficient of damping rate for blade lag damper. 

thrust coefficient of the rotor; thrust nondimensionalized by 

(p7rR^2). 

plate bending stiffness. 

Xg coordinate of coincident flat-lag hinge or hingeless blade 
offset point. 

chordwise position of section tension center, (+) forward from 
elastic axis. 
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LIST OF SYMBOLS (Cont’d.) 


Elyj Elg 

[e], [pj 
fa 

{p(«/»)l 

^PR 
F • F 

y?’ ^5 

g 

GJ 

W 

^ 10 » ^^10 


%F 


element of E and P matrices, respectively. 

blade bending stiffness for flatwise and edgewise bending, 
respectively. 

torsioneO. transfer matrices resulting from cascade multipli- 
cation and other algebraic manipulations . 

spanwise angle of attack correction factor to accovint for 
three-dimensional tip effects. 

excitation vector in time-history solution fozm of nonlinear 
dynamic equations. 

push- rod load, (+) when directed in the direction. 

concentrated forces applied to the blade in the y^ and 
directions, respectively. 

acceleration due to gravity. 

torsional stiffness of blade section. 

identity matrix of dimension m. 

radius of gyration of section tension carrying area. 

mass radii of gyration of blade section about axes through 
and perpendicular to the spanwise (x^) axis and in the , 
chordwise and thicknesswise directions, respectively. 

spring rates of edgewise and flatwise retention springs, 
respectively, of torque tube to blade spar at flexbeam - torque 
tube juncture. 

induced velocity gradient factor. 

spring rates of torque tube principal axes (edgewise and 
flatwise) retention springs, respectively, of torque tube to 
flexbeam at snvdiber end of torque tube. 
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LIST OP SYMBOLS (Cont'd.) 




K01 


spring rates of inplane, and vertical retention springs, 
respectively, of torque tube to flexbeam at snubber end of 
torque tube. 

spring rate of flexbeam torsion due to elastic and bifilar 
effects . 


K02 


spring rate of effective root torsion restraint due to control 
system flexibility and flexbeam bending. 



spring rate of root torsion restraint due to control system 
flexibility. 


i 




m 


length of element of flexbeam over which torsion properties are 
assumed constant nondimensionalized by flexbeam length. 

length of beam element over which the tension is assumed 
constant. 

vector of shears, moments, and deflections given as linear 
combinations of other specified deflections. 

subvectors of | , having shear, moment, and deflection 
components each in y- and z- directions, respectively. 

blade mass distribution. 


m^ 


reference blade mass distribution, taken to be that of the 5th 
blade segment. 


M 


flatwise bending moment at an end of a beam element; also Mach 
number. 


Mu) 


moment due to lag damper . 

blade root torsion restraint moment. 


*V5» ^5 


components of blade moment about axes in the 5rCoordinate 
system, arising from inertia, aerodynamic and/or concentrated 
mechanical effects. 


HEM, NPM, NTM numbers of assumed flatwise, edgewise, and torsion natural 
"uncoupled" primitive modes, respectively. 
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LIST OF SYMBOLS (Coiit’d.) 


HH mmiber of haxmonlcs in vajrlable inflow description. 

HSEG number of segments into idilch the blade Is divided each having 

assuTOd spantilse ccmstant section properties. 




section shear load distributions in directions of axes in the 
coordinate system. 


<bc 5 » Qyjj <3355 
LQjtJ 


general e3^>resslon for a response variable deflection, 
element of vector. 

blade k'th edgewise modal response variable. 

blade i'th flatwise modal response varlcd>le. 

section moment load distributions about axes in the 
coordinate system. 

blade j'th torsion modal response variable. 

vector of torque tube deflection effectivlty constants. 

blade spanwlse coordinate, measured from offset, e, in : 
direction. 


R 

•ij 

S 


rotor radius. 

Intermediate torque tube bending stiffness matrices, 
element of the Sq matrix. 

flexbeam length (span), or flatwise shear at an end of a 
beam element, as appropriate. 





caapaaenta of concentrated shear in directions of axes in tlie 
coordinate system. 



trigonometrically' resolved torque tube bending stiffness 
matrices . 




element of a T matrix. 





LIST OF SXMBOIS (Cont'd.) 


T tension at an arbitrary blade spanwlse station. 

Tji^ constants for polynomial representation of spanwise variable 

torque applied to flexbeam . 


l\] 

[■'ll 

[t5] 


^ elastic transfer matrix at k*th semi -segment. 

{^l]> [^3] »[^]» coordinate system transformation matrices relating rotating 
coordinate system deflections to the inertial frame. 


"e 

u 



elastic spanwise deflection of arbitrary spanwise coordinate. 

resultant air velocity relative to blade section, 

air velocity component relative to blade section in (+) 
direction. 




U 

X 


air velocity component relative to blade section in (-) 
direction. 


air velocity component relative to blade section in (+) x^ 
direction. 


Vg, Wg elastic edgewise and flatwise bending deflections, respectively, 

of an arbitrary spanwise coordinate. 


Vj^ , Vi , Vi zeroth and n'th cosine and sine components of variable induced 
o Ug Ug velocities . 


Vp, Vj 


components of air velocities relative to blade section devoid 
of cosine sind/or sine of total pitch angle and comprised of 
response variables. 

components of air velocities relative to a blade section in 
chordwise and thicknesswise directions, respectively. 


Av, Aw 


deflection correction terms due to first order twist effects. 


V strain energy. 

AV, AW deflection correction terms due to second order twist effects. 
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LIST of SYMBOIS (Cont'd.) 


Xg, 72, Z2 
X5, 75, Z5 
Ax 

{* 5 } 

7io» Zxo 

A A 

yiO> ^10 

yi0e/4> 

^lOCO 
{ay}, {az} 

Az 

z 

“eff 

aqs 

“r 

)8 


components of the 2-coordinate S7Stem, defined to be affixed to 
the rotating hub, 

components of the 5-coordinate S7stem, defined to be rotating 
with the hub, but at the blade coned and lagged position. 

length of blade spanwise segment over which the spanwise 
properties are assumed to be constant. / 

components of blade deflection in the inertial and 5-coordinaie 
S7stems, respectivel7, as measured in the 5 -coordinate' S7stem. 

chordwlse and thicknesswise position coordinate, respect! vel7, 
of an arbltrar7 point within a blade section. 

Inplane and out-of -plane position coordinates, respect! vel7, of 
an arbitrar7 point within a blade section. 

chordwise distances of quarter chord, three quarter chord, and 
mass center, respectively, from the reference axis at a blade 
section. 

subvectors of |a^ , having specified deflections in 7 and z 
directions, respect! vel7. 

deflection in z^ direction of the push-rod attachment point 
due to an7 and all blade bending deflections, flapping and 
lead-lagging, with the push- rod disconnected. 

flatwise (vertical) deflection per \mlt force of inboard end 
of cantilevered torque tube, due to flexbeam bending 

effective blade section angle of attack approximately corrected 
for three-dimensional tip effects and used for evaluation of 
section aerodynamic coefficients. 

quasi-static blade section angle of attack. 

rotor (hub axis) angle of attack. 

blade flapping (or precone) angle; also used to denote tension 
(or compression) to elastic bending stiffness parameter for 
a beam element. 
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'>'0RB 

rBPej 

8 

{4 


V 

A 

■^k, k+1 


deflection mode shape for the k'th edgewise normal mode, 
deflection mode shape for the I'th flatwise normal mode, 
deflection mode shape for the j'th torsion normal mode, 
deflection shape for the pseudo- (rigid body) torsion mode, 
blfllar-twist coupling function for j'th torsion mode, 
blade lead angle; also used to denote perturbatlonal quantity, 
vector of specified deflections. 

small nuoiber (less than unity) used to assess relative orders 
of magnitude for various terms. 

flexbeam spanwise coordinate, measured from offset, 
nondimenslonalized by S. 

spanwise coordinate of boundary point between k'th and k-fl'th 
segments . 

built-in blade section twist angle; i.e., that section pitch 
angle resulting when the aerodynamic pitch angle at 75% span . 
is zero deg. 



bvdlt-ln twist angle of flexbeam measured when flexbeam is 
torslonally unloaded. 


6, 6q elastic torsion deflection angle. 

blade root torsion deflection angle. 




abbreviated torsion flexlbllitv coefficients containing 
deflection . dependency. 


Oiji 


blade tip torsion deflection angle. 


^vx» automatic blade pitch change per unit deflection In k'th 

6^ edgewise mode, I'th flatwise mode, flapping and lead- lagging 

motions, respectively. 
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blade pitch dngle due to input control angle. 

Ad automatic pitch changes accruing from any and all blade 

bending deflection, flapping and lead- lagging. 


e 


total local blade pitch angle. 


A0J 



flexbeam torsion flexibility coefficients. 

total torsion deflection of flexbeam at jvinction point beyond 
built-in vsQue. 

blade section pitch damping eff activity factor. 


X 


total inflow ratio, or section torsional characteristic constant, 
as appropriate. 


uniform portion of Inflow. 


^RAM 


inflow contribution from forward flight "ram” effects. 




advsuice ratio. 


V 


Poisson's ratio. 


p air density. 

(f> inflow angle . 


i/r blade azimuthal (angular) position. 

23 (nondimenslonal) frequency used in queidrature formulae, talcen 

to be \incoupled nattiral frequency of degree of freedom. 




(nondlmensional) uncoupled natureO. frequencies of i'th f^latwise 
bending mode, k'th edgewise bending mode and J 'th torsion mode, 
re spec tively . 


ft 


rotor rotationail frequency. 
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LIST OF SYMBOLS (Cont’d.) 


Subscripts 

( )a 

( h 

< )d 

( >. 

( )e 

( >F 

( )pB 

( )g 

( )l 
( )j 

( )o 

^ w 

( )s 
( 

Superscripts 
(~ ) 

( ) 

( )’ 


effects of aerodynamic origin. 

structurally built-in pea-ameter, or conditions of blade 
immediately outboard of Juncture. 

effects of dynamic origin. 

due to elastic deformation. 

in edgewise (section major principal aixis) direction, 
in flatwise (section minor principal axis) direction, 
flexbeam . 

effects of gravitational origin, 
conditions at inboard end. 

conditions at flexbeam -torque tube Juncture, 
conditions at outboard end. 
push- rod. 

conditions at the inboard end of the torsionally active 
portion of the blade. 

conditions at snubber (inboard end of torque tube), 
torque tube. 

nondimensionalization by combinations of m^, R and/or £2. 
differentiation with respect to (ftt). 
differentiation with respect to (r/R). 
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LIST 6F symbols (Cont’d.) 



pertains to loads arising from deflections . 


)(ER) 


pertains to loads arising directly from push- rod load. 

pertains to first and second parts of definitions of 
deflection correction teams, respectively. 


)H( )W 


to be evaluated at a specified point minus or plus an 
infinitesimal amount. 
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PRINCIPAL ASSUMPTIONS 


The principal assumptions used to derive the basic differential equation 

of motion are as follows: 

1. The rotor is rotating at a constant angular velocity, has infinite hub 
impedance, and is in steady translational flight. 

2. The blade elasticity is adequately described by| the conventional (linear) 
beam bending and bar torsion characteristics described in Ref. 4. Al- 
though the effects of the additional section constants and B 2 described 
therein are usually considered to be negligible for helicopter applica- 
tions, they are potentially important for accurately analyzing solid 
sectional, highly twisted propeller blades and/or wind turbines. To pre- 
serve consistency with the rigor applied to other aspects of structural 
twist and to achieve universality with such nonhelicopter rotor systems, 
these elastic section constants are retained in the full nonlinear formu- 
lation given in Ref. 4. 

3 . The elastic (torsion) axis of the undeflected blade is a straight line. 

•However, when deflected in bending, the elastic axis defines a space 

curve about which the local torsion deflections must take place. 

4. The blade aerodynamic and structural twist distributions are nonlinear; 
additionally the structural twist of the flexbeam (bearingless rotor 
applications only) is time variable. 

The total (Integrated) angle of structural twist is negligible beyond 
second order; cases of large local twist rates over short sections of 
span are not dented, however. See section on coordinate system for more 
details . 

6. Radial foreshortening of blade elements due solely to elastic deflections, 
in the absence of precone (or flapping), and prelag (or lagging) is ade- 
quately represented by a second order function of flatwise bending. 

7 . The feathering axis is coincident with the elastic axis of the elastically 
undeformed blade . 

8. The blade distributions of center of gravity, aerodynamic center and 
center of tension (intersection of flatwise and edgewise neutral axes) 
are, in general, noncoincident with the elastic socis. 


9> The blade sections have finite thicknesswlse mass, but the thlcknesswise 
displacements of the section center-of -gravity away from the chordwlse 
principal axis is negligible. 

10. Vfhlle assumptions regarding the smallness of various quantities and 

products of these quantities are not generally required for the ioqple- 
mentation of time-history solutions of the full nonlinear equations, they 
are required for effecting consistent linearized approximations for the 
elgensolutions. For this case, coefficients of the perturbatlonal varl* 
ables, whose orders of magnitude exceed are neglected. Here e is an 
unspecified small number less than unity and where the assumed orders of 
magnitude of the various pertinent quantities, as measured by e, are 
given in Table I . 
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COORDINATE TRANSFORMATIONS DUE TO PITCH ANGLE AND IWIST 


The aeroelastlc analysis Is a modal type analysis 'which uses, as Its 
basic description of the elastic deflections, hereafter called primitive modes, 
the normal vibration modes calculated assuming the blade to have zero coning, 
pitch angle and twist. These modes are sometimes referred to as "uncoupled” 
modes. Thus, the resulting aeroelastlc bending responses, as ultimately 
coupled by coning, pitch angle, twist, aerodynamics, etc., must be Intezpreted 
as flatwise and edgewise responses rather than out-of-plane and Inplane re- 
sponses. The basic advantage claimed for this choice of primitive modes over 
those appropriate to a blade already twisted and pitched at some nominal angle, 
which are sometimes referred to as "coupled" modes. Is the convenience of usage 
especially for bearingless rotor applications. A single set of uncoupled 
primitive modes needs to be calculated for any one rotor speed and Is adequate 
for all subsequent twist and pitch angle combinations. The purpose of this 
section, therefore, is to describe the introduction of pitch angle and twist 
using uncoupled normal modes as the set of primitive modes. 

The basic blade coordinate system used herein and referred to as the "5" 
coordinate system consists of the elastically undeflected blade at the flapped 
(or preconed), lagged (or prelagged) position. See Ref. 1 for a detailed dis- 
cussion of the coordinate transformations preceding this. As shown below In 
Fig. 2, the x^-axis is out the span from the blade root or offset point, as 
appropriate; the yc-axis is parallel to the hub plane and positive forwai^i, 
while the z^-axis Is perpendicular to x^ and y^ and positive upward, but not 
generally parallel with the spin axis. 

As Is discussed above, the local elastic deformations (nondlmensionalized 
by rotor radius, R) consist of finite series summations of normal bending 
modes In the flatwise and edgewise directions: 


NFM 

i = l I I 


NEM 

V. = E r m m 

k=i « ^ 


( 1 ) 


( 2 ) 
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Figure 2, - Schematic of the ''5" Coordinate System. 


In the presence of only blade pitch angle the resulting deflections in 
the "5” coordinate system are simple trigonometric transformations of these 
flatwise and edgewise contributions. The effect of twist, however, is to 
require an "Integrated" trigonometric transformation. This integrated effect 
can be achieved by means of , first , a direct trigonometric transformation not 
on deflection, but upon second (nondimensional) spanwise derivatives of the 
deflections : 


Vj s Vg cos® - Wg sin® 


(3) 


Zj = Vg sin ® + Wg cos® 
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This coordinate system transformation has the advantage that the "force” 
boundary conditions at the tip of the blade in the "3" coordinate system are 
always satisfied. Equations (3) and (4) and their first spanwlse derivatives 
taken together with the boundary conditions imposed on and Yvv given in 
Eqs. (l) and (2) axe sufficient for this result. Eqs. (3j and ( 4 ; can then be 
integrated by parts to give the fundamental deflection coordinate transforma- 
tion used throughout the g 400 aeroelastlc analysis. This transformation 
becomes the usual trigonometric transformation on deflections given in Refs. 1, 
4 and elsewhere in the literature, but with the addition of various deflection 
correction terms due to twist: 


Yg = (7g+Av-AV)cos@- (Wg-Aw-AW)sin© + 0(®[^ (5) 


Zg = (Vg + Av -AV) sin© + (Wg-Aw-AW)cos© +0(®^ (6) 


where. the underlined terms are, by assumption, negligible, and where the 
deflection correction terms are given by the following: 

first order in twist: 


Av = [ di; 


= (Av['^+ Avj^^q^. = AVjQ^. 


(7) 




'o'o 


( 8 ) 


= (Aw"’+Aw,;'’)q =Aw^q 
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second order in twist: 


= [/e'Aw„dr, + rV 0'^wf drjdfJ q. 


( 9 ) 


= (Av‘'> + ^v'®)q = AV^q 


'O'O 

k 



= (Aw['' + Awf^)q^^. = AW.q„, 


It is to be noted that the total ( nondimens ional) twist rate, 0' , contains 
the built-in twist, 6g , the twist due to control inputs, 0^ , and the time 
dependent elastic deflection, 0 g Thus, the deflection correc- 

tion terms, Av, Aw, AV and AW, nominally contain nonlinear products of q^j 
with either and/or Herein, however, the products qg. • q^, and 

q0. • qy^ are retained only in the Av and Aw (first order) correction terms and' ' 
denoted as Av^ and Awg respectively; the AV and AW (second order) correction 
terms retain only the built-in twist rate, 0 g , and that due to control Inputs, 
0 , in accordance with the assumed relative order-of -magnitudes given in 

Table I. 

TO complete this section, a physical interpretation of the above formulated 
coordinate transformation is presented by considering the following argument: 
Let a uniformly twisted blade possessing only flatwise flexibility be uniformly 
bent (in only the flatwise sense) over its span by continuously taking small 
bends at points outward along the span starting at the root, as shown in 
Fig. 3 below. Then, with each such small bend, the deflection locus of the 
blade tip, point P, is traced to its final deflected position, point p’ . The 
initial part of the locus must be in a direction normal to the blade root 
chord since the initial flatwise bend is defined to be normal to the blade 
root chord. Similarly, the final part of the locus must be in a direction 
normal to the blade tip chord since the final flatwise bend is defined normal 
to the blade tip chord. Thus, the locus must define a curved path (arc PP'), 
as shown. In contrast, using the straight trigonometric transformation on the 
flatwise deflection, Wg, together with the pitch angle at the tip, 0^, the 
blade tip would be predicted to be at point P", where the flatwise deflection 
equals both the PP" and the arc length PP\ The figure clearly shows that the 
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Figure 3‘- Illustration of the Physical Significance of the Deflection 
Correction Terns Resulting from Assumed Coordinate Transformation. 


straight trigonometric transformation will yield the actual deflection only if 
the chordwise and thicknesswise corrections, QP and QP", respectively, are 
added to the trigonometric transformation; 



Yj = ( 0 P^) COS - { Pf>"- 0 p") si n dj 


Zg H (Qp') sin By + ( PP'- QP") cos^y 


which reduce to Eqs. (5) and (6) where PP" is given by Wg and where QP* and 
QP" are, respectively, approximated by Av and AW. 
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DERIVATION OF BASIC CmMIC EQUATIONS OF MOTION 


Flatwise and Edgewise Bending Equations 

Hie differential equations of blade bending are obtained by equilibrating 
the flatwise and edgewise moments at arbitrary span locations. Reference to 
Fig. 2 and the principal assumptions leads to the following moment equilibrium 
equations; 


Elyw;' = -MygCOsQ-MzjjSin® (11) 

ElzVg- e^T -EBg = MjgCos® -Myosin© (12) 


. Application of the Galerkin technique requires that the equilibrium 
equations be in the form of loading equations; hence, Eqs. (11) and (12) must 
first be doubly differentiated. A method of performing this double differen- 
tiation which maintains equation compactness is given below. The first differ- 
entiation yields the following equations; 

flatwise bending ; 


[ElyWe] +® [ EIjv;'- e^T -EBg ( 0 B+I^e)^e] = -[w'ygCos® + M^^^sin®] (13) 


edgevise bending ; 

[ei^v" - e^T -EBg (fi^+^5e)®e ] ' = [m'^^cos® - Myosin®] (l4) 
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where the derivatives of the applied moments , mI and are expressible as: 

^5 5 

My^=^ pj^dri-ZgT-qy^ (15) 


”' 5 ' + - ( 16 ) 

and where the tension Is given by: 

T=|\dr, (17) 


The applied force loading distributions, and and the applied moment 

loading distributions qy_ and q 2 _, all Include not only the Inertia and aero- 
dynamic contributions bu^ any loads due to mechanical loads applied by push- 
rods, dampers, etc. Ihe second differentiation of the moment equilibrium 
equations Is performed after first resolving Eqs. (13) and (ih) to Isolate 
1^5 and 


-My^= { [(EIjV^- e/,- EBg ®EIyW;] Sine 

+ [(Eiyw;)+e'(Eijv^- e^T - EB2 (04 ©;)] cos®}' (rhs)^ 


(18) 


{[(EIzV^e^-EBa (©' +^0')e;)'-®'EIyw;]cos® 
-[(EIyW;)'+®'(El2Ve-e/-EB2 (fi^+^0e)®e^]sin®}' = (RHS)^ 


(19) 


Application of the Galerkln technique requires consideration of the work 
due to virtual displacements of the bending modes. It must be noted, however, 
as was shown In the preceding subsection, that In the presence of twist, flat- 
wise and edgewise modal displacements each generate deflections In both the 
flatwise and edgewise directions. Hence, the proper formal application of the 
Galerkln technique Is to Integrate the Inner products of loadings and virtual 
displacements : 
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(20) 




k*th edgewise modal equation ! 


(31) 


The quantities (RHS) and (RHS)^, as defined by Eqs. (l8) and (19) > involve 
differentiated conib inat ions of twist and cosine and sine of the pitch angle 
and elastic restoring moment. Although not immediately obvious, the integrals 
of these (RHS) quantities with the components of deflection, as given by 
Eqs. (5) and (6), can be evaluated, using integration by parts, to yield the 
following siu 5 )lified forms; 


R R 

I yjEIyW;'dr (22) 

R R 1 

/ [(z,) (RHSj^+Cy^) (RHS(,]dr = / (23) 

OR R Or 


The desired basic bending modal excitation equations can then be written by 
combining the results of Eqs. (5), (6), (15), (l6), (20), (21), (22) and (23). 
Hereafter all quantities and equations will be written in nondimens ional form 
(without overbars) where the nondimens ionalizat ion is accomplished using appro- 
priate combinations of R, fi and m^: 
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cose-pt, sin®)-Av,(pn sin® + pu„ cos®) 

■b » ^5 Vs ' ^5 V5 

+ T[(r^- AwJ^^Kw^- + Av/^^(v4+ 

+(y^-Aw[®Kqoy^cos® + Qo^^sin®)- Av[®(qo^cos®-qoy^sln®) + y^EIyWe}dr 
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The subscripts A and D denote terms of aerodynamic and dynamic origin, 
respectively, and the finite summations over m represent the modal excitations 
due to a finite number of concentrated forces and moments. 


Torsion Equation 

In a development similar to that given above, the differential equation 
for torsion Is obtained by first equilibrating the blade torsional moment. In 
Ref. 5 Is derived the torsional equilibrium equation for the blade, consistent 
vlth the assumption of a space curve torsional axis: 


GJe;+0'k2T + ^EB, e'-EB^'eive 

= /{[‘^-^5)y^(y5,-y5^z'5]Px5(r,) 

+[(y5,-y5)-^-»')y5] Pzs^n) +qxs(n^ +ysqy (r,) +2^q (r,)}dri 


(26) 


where y^^ and 


are the Inplane and out-of-plane deflections evaluated at r^^. 


Differentiation of this equation yields an Intermediate form of the 
required torsion loading equation: 


[gj 8;+ efkjT+-^EB, •' -EBjIfliv;']' -- y^q,^- z^q,^ 

■ 2?jT' [(>5, - I's* Pzjl^ > - (f| -''Pyjlfl ' - Iz, (n )] <*1 
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Integration by parts yields a second intermediate form; 


[6J0;+@'k|T +^EB, = -q,^- 

+ ^5 jf ks, jT drg -jf^ Pzg(r2) drg + qy^(r, )] dr, 

" ^5 / bs, / Pxg^'z )dr2 -/ p (r2)dr2 - qj (r, )] dr, 

n *"i r, ^ 5 -• 


(28) 


which, when coiribined with Eqs. (15), (l6), and (17)j with (11) and (12), and 
finally with (3) and (4) yields the final desired form of the loading equation: 


[g JO' + ©kjl + EB, (e'^ - e'- EBg^; v" ] 

- [{ EIj- Ely ) vX- ( e^T + EBg (e' + -^^e) ) Wg] 

Application of the Galerkin technique to this equation then follows in a 
straightforward manner to give the following basic form of the J 'th torsion 
modal excitation equation: 


(29) 


= “dxg 


hox." ysPov" ^sQdz + (El2-EIy)v^w4'-(e/+EB2(e;+^0'e)5;)w^] 
0 J ^ '5 5 

+ Ye. [gj 9; + O'kjT + ^ EB, (e'^ - e^) 9' - ]} dr 



d'’ +E’'d. ('m)M 

5 m=l J 



(30) 


where, again, the summation over m represents the effect of concentrated torsion 
loads on the torsion modal excitation. 
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Rigid Flapping and Lagging Equations 

Ihe equations of motion for rigid flapping and lagging are obtained by 
equilibrating the moments about the articulation hinge about the and 
axes, respectively; 

rigid flapping ! 

- 0 (3D 

rigid lagging; 


ro-Pys-VsPxs+Qz,'-!' 


+ “uo 


m= I 



= 0 


(32) 


•where is the root moment due to the lag damper and is typically proportional 
to blade root angular rate; 


M 


LD 



(33) 


with Ctjj as the familiar linear damping rate of the damper. Note that as per 
Eq. (5)the ( nondimens ional) time derivative of elastic in-plane slope will 
generally contain contributions proportional to flatwise and edgewise modal 
deflection and to flatwise, edgewise and torsion modal rates. 

The equations of motion presented above for the bending and torsion modes 
and for rigid flapping and lagging ((24), (25 ), (30 ), (31) j and (32), respec- 
tively) are complete only in so far as the load distributions are (explicitly 
or implicitly) available. The following section presents implicit statements 
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of the loading in terms of the blade deflections, as formulated in the 
following section. Furthermore, these equations are basic required forms of 
the aeroelastic response excita*bion equations from which, together with appro- 
priate expressions for the load distributions, a set of linearized equations 
can be expanded for eigensolutions, or a set of nonlinear equations of the 
form; 




( 34 ) 


can be written for time-history solutions. The A matrix above consists of the 
coefficients of the second derivatives of the modal responses and represent 
terms extracted from the compact expressions for (Synamic load distribution. 
This matrix is also discussed more fully in the following section. 
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DYNAMIC AND AERODYNAMIC LOAD DISIRIEOTIONS 


The load distribution a;^ear^g in the above derive dynamic equations for 
the response variables were expressed only implicitly for two reasons: first, 

because it simplifies the explanation of the application of the Galerkln tech- 
nique and second, because ccmpletely expanded, explicit, expressions for the 
loadings are tedious, cumbersome and not actually required for the more impor- 
tant time-history solution. Complete linearized expansions of the loadings 
are required, however, for a formulation of the linearized equations used in 
the elgensolution. 


Dynamic Load Distributions 

Theusuai approach of using D'Alembert 's principle to express the inertial 
acceleration of the distributed blade mass as an equivalently applied dynamic 
load distribution is followed herein. The position vector of a point mass 
particle, with a chordwise and thicknesswise displacement relative to the "5" 
coordinate system, is written as follows: 


I r+Ug-yj^^VgCos© + ^sin©)-Z|Q(^cos©-y5sin©) 
{^ 5 } = < y5 + ijoCos©-z^sin© 

V ^5 +z,qCos© 


where y^^Q and z^q are, respectively, the (forward) chordwise and (upward) 
thicknesswise locations of the point mass from the reference, x^, sixis, and 
where the axial deflection due to elastic flatwise bending is given by the 
following expression: 



The displacement of the point mass particle relative to the inertial 
frame is obtained by means of four consecutive coordinate transformations and 
is written in the following compact form: 
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{x,} = [t,] 


+ T, 


JW{x,}] 


(37) 


where the four coordinate transfoxmations, described in detail in Hef . 1, in 
consecutive order, account for rotor rotation, blade root offset, lead-lag 
rotation about the (articulation) hinge or offset point, and flapwise rota- 
tion also about' the offset point: 


fcos^ 

-sln^ o' 

[T,J = sln^ 

cos^ 0 

1 0 

0 1 . 


(38a) 



cosS 

sinS 

0 


-sinS 

cosS 

0 


0 

0 

I 


(38b) 


cos/3 0 -sin/91 (38c) 

0 ‘ ® I 

sin/3 0 C0S/3J 

Upon letting: 





(384) 


the inertial acceleration of the point blade element can be written as: 




( 39 ) 
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^idiere the details of fonalng the indicated differentiations and matrix 
multiplications are omitted herein for clarity. The dynamic force load dis- 
tributions are formed as follows: 

S) 

So, > (to) 

q ® \ blode 
/ section 
® area 



whose components are formed by using Eq. 


(39) aJid can be written as: 


Pdxj" nn{e+r(l+2§-/9^) + 2jf5 + Ue 


(41a) 


pDyg= -m{eS + r(S-2/34) +y5+2(Ue-)0^-425) -(l + 28)y5 


+y,0cg[*2( *e - + Av^^’) - ( 1 + 2§ + ®^+2/3|» cos© - (Si + 24)sin©]} (4lb ) 


■f^{®^‘*T(^*+)8+2/38) + Z 5 + 2 i 8 y 5 + y,Q^^[®cos®- (^^+2/30) sin®]} (4lc) 


The dynamic moment load distributions are similarly formed: 



-ffpdA 

blO(je 

section 

area 


0 i|o 

-Zio 0 
Vio 0 



0 


(42) 


where the chordwise and thicknesswise integration variables resolved into the 
”5" coordinate systan are given by: 



^io=yioCos0-ZKiSine 


(43a) 


^K) = Viosind+Ziocos© 

The conrponeiibs Qf the moment Xoad distributions can then be written as 

^0x5= + 

-[e8 +r(s‘-2/8/9) +‘^-2(*e +)9?5+i3z5)-(l + 28)yJsln®} 

+ ky^^{® + 2 /Scos^@ -(i+2i-/S^)sin®cos0 

+ 2( 'Xg - cos®} 

+ kf |® + 2/Ss|n®® + (i+28-;3®)sin®cos® 

*10 L 

+ 2{v 4 - AV^^'+Av^^^)sin®}^ 

qDyj^=-"'^-y,OcG{e+^('+2S-^^) + 255+Ue-/3z5}sin® 

+ k 5 io{(l+ 28 ) [v/g- AW^^^ - Aw^ V/3C0S®] 

+ 2§[(l + §)cos® +4s>n®]}cos® 

+ kf|o{(l+2i)[v4-Av‘^^+Av^'+)9sin®] 

+ 2®[(I+S) sin®-4cos®]}sin(^ 


(43b) 


(44a) 


(44b) 
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+ kyio{( l+2§) [wg- Aw‘^^V/8cos©] 

+ 2 <§[(l +f )cos® +)§sin®]]- sin© 

- k||o{(l+2S)[y^- Av^^^+iSsin®] 

+ 2 ® [(I +§) sin®-/8cos®]} cos(S^ 


where ^10 respectively, the chordwise mass center and 

the thlcMesswlse and chorawise mass radii of gyration. 

For both the eigensolution and the time history solution, those dynamic 
load terns containing second time derivatives must be extracted to form the 
inertia coupling matrix, [A], indicated in equation (3^). Using the load 
distributions given above, the elements of this matrix can be conveniently 
written using the following partitioned representation of the A matrix: 






Awj/3 

AwjS 

Ay^Wj 


\»i 


\s 

Vi 

Vk 



AdjS 

A/5W. 

A^v^ 

V 


0 

AfiWj 



0 

AfiS 


where: 


^WjWjn” f +AVj AVp^jdr 


(^5) 


(46a) 


AwiVjj = Av^Wj 


J‘m[(yy -AV,^) Av. - (r^- AWj)AwJdr 


(46b) 
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T1 '0 


NF»* 





<’A 






K'--' 


^y<.^ * '0 

^0 


. r'mr\iy^-^'^-. 

•*0 

.. rU\.^V.,CDS® 

^r\ 




-A*,!*"®!'*' 


-tv'""! 


(U6d) 

(U6e) 

(U6f') 




> 


»v^»„'Jo 

/; •rt 


N* i '0 






IV - ^V^^^^'**ll/ 


vtWV"'‘'^'“‘''’' 

v'»\ 






(U6g) 


(U6la) 


A,'p - '^'‘>1 -0 


•ri 


C0S®V^ 


(U6i) 


^v^^S * '0 


.r'm-lv 













- AWj)dr 








(46k) 




r)^ cos® dr 


(461) 



(46m) . 


(46n) 


Aerodynamic Load Distributions 

The aerodynamic load distributions used in the analysis are assinned to be 
two dimensional and the usual strip-theory techniques as typically described 
in Ref. 1 are therefore employed. More specifically, at each spanwise station 
the two-dimensional airfoil section angle-of -attack is calculated based upon 
section geometric pitch angle, 0 , and inflow angle, 0, based upon airflow 

velocities at the 3/4 chord point? 


Qqs= @ + 4 > 


® + tan”’ 



(47) 
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where ; 


^ = x(l--^) -U^[zj(l--^)+/3j-!5-r4-y,o^(5+®'Ux)cos® (48) 

4 

= e +/xsin('/'+8) + r (i + Uxyg^i--^j 

+ yg + rS-y,Q {5 + ©'Uj() sin® 

T 

where the radial velocity component is given by: 


= ^cos{'/'+8) + y,o^^cos® (50) 

T 

and where the derivatives of the "5" coordinate elastic deflections are obtain- 
ed from appropriate differentiation of Eqs. (5) and (6). The total rotor 
inflow, x(r,i|/), is comprised of a "ram" portion due to forward flight and a 
(harmonic) induced velocity portion which is assumed to be,, in general, both 
radially and azimuthally variable: 


« = ''RAM - -E [»i Jr)cosnt + v.^^(r)sinnV'] 


(51) 


where : ■ 


V 

^RAM = sinap = /xtanOp 


(52) 


Equation (51) represents a completely general description of the harmonic 
rotor inflow. Regardless of how the induced velocity components are obtained, 
either experimentally or analytically, they are accepted by the analysis as 
an environmental excitation and used directly as per Eq. (48). 

Two specific, optional assumptions which can be made on the induced 
velocity are as follows: 


39 



a. 


unifoim inflow: 


V, (r) = — . • 

'0 


where \ 


n = l, 2 , 


ns . 


b . generalized Glauert (momentum) : 


V. (r) = (uniform) 
lo o 




V. (r) = rv,, 
'is ' 


V; (r) = V. (r) = O 5 n =2, 3, , NH 


■nc 


*ns 


^0 ' ^RAm“ ''0 


(53a) 

(53b) 

(54a) 

(54b) 

(54c) 

(54d) 

(54e) 


The nonzero (zeroth and first harmonic) components of induced velocity can 
then be related to rotor steady thrust and hub moments using momentiun consi- 
derations : 



(55a) 

(55b) 


(A^^+4-^oVo) (55c) 

The Glauert induced velocity gradient factor, K^, in Eq.. (54b) ‘is used to 
account for the strong cosine ccmponent of inflow present even for zero pitching 



moment. The approximation to used herein Is that given in Ref. 6 based 
upon results of Ref. 7: 


Kv = 


4/1 

3.6 |Xol+3/x 


( 56 ) 


The foimulations given above for the local velocities together with Eq,. 
( 4 ?) is sufficient to define a quasi-static two-dimensional angle-of -attack. 
Approximate three-dimensional tip effects are introduced by loultiplying this 
quasi-static angle-of -attack by a function 'which is unity over most of the 
blade> but reduces abruptly to zero at the blade tip: 



0 < r < 2B-I 


2B-I < r < I 


(57a) 


where B is the conventional tip loss factor. 


The effective angle-of -attack used to calculate aerodynamic coefficients 
is 'then given by: . 


*eff 


ta®qs 


^aaqs'^®9n(aqs){l-fa)^$ 


‘qs 

*qs 


^ -f 


lOcsl > -Y 


(57b) 


The Mach number at the airfoil section is given simply by the following 
e:qoression : 


M 



( 58 ) 


From the effective angle-of-attack and Mach number, the following e:g>res- 
sions for (nondimensional) aerodynamic load distributions can be formed: 


4l 



(59a) 


(■^'o') ^d'Jp> (59b) 

"‘Xs" H/ 4 ('’Az 5 “S®-% 5 Sin®)+ i CK.^S] (59c) 


where : 


II 

0> 



( 60 a) 

Cd = Cd(aeff'^) 



( 60 b) 

^ ^nnc/4^“eff> 



( 60 c) 


lOqsl 

lotqsl 

>1- 

( 60 d) 


ACjj = incremental drag coefficient introduced to account for surface rough- 
ness. 


The dynamic equations given in the previous subsection together with the 
load distributions presented above are sufficient to complete the basic aero- 
elastic analysis of the rotor blade. Such a basic analysis, however, omits 
the effects of pitch-f lat/edge coupling. The following subsection incliades a 
unified method for including such coupling. 
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RIGID BODY FEATHERING MOTION AND 


BASIC PITCH-FLAT/EDGE COUPLING EFFECTS 


The simulations of* the rigid body feathering motion (if present) and the 
automatic blade pitch changes accruing from elastic bending deflections and/or 
flapping and lagging motions are both acconrplished by the introduction of a 
torsional "pseudo "-mode. Such a mode is no more than the rigid body feather- 
ing motion of the blade as would be generated by a control input. Operation- 
ally, this "pseudo"-mode is treated in the analysis as a spanwise variable 
mode shape in addition to the conventional normal torsion modes used to 
describe the blade torsion. In general, this mode has a unit value over the 
blade span except, for analyses of the bearingless rotor, over the flexbeam - 
torque tube portion of the blade wherein the spanwise elastic torsion deflec- 
tion of the flexbeam (due to a control input) is used. Reference to Fig. 4 
shows the pertinent feat\rres of this pseudo-mode, especially as it is applied 
to bearingless rotors: 



nonoimensional span , f 

Figure 4 . - Details of (Rigid Body) Torsion Pseudo Mcxie. 
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Note that the elastically effective portion of the pseudo-mode can, 
optionally, be taken to be a linear function typical of rod- like torsion or 
an ogee function more typical of plate-like torsion. Finally, it should be 
stressed that this pseudo-mode is not, in general, orthogonal to the blade 
normal modes in torsion. 

For bearingless rotor applications, the pseudo- torsion mode serves the 
unique function of providing a convenient vay of introducing time-variable 
structural twist (l.e., due to control inputs) into the analysis. Ey treating 
this component of twist as but an additional torsion mode proportional to con- 
trol input, the following expression for total pitch angle can be written; 


NTM 


0 


j'l j j 


Similarly, the total twist rate can be written as; 


■NTM 


(6la) 


(6lb) 


Sy using this expression for twist rate, together with the appropriate 
pseudo-mode, the Incremental deflection functions, Av, Aw, AV and AW described 
in an above subsection, can be completely formulated. 

Ihe pseudo-torsion mode serves two additional functions common to all 
types of rotors: First, it provides a convenient way to include a rigid-body 

torsion degree-of- freedom arising from root torsion restraint (control system) 
flexibility and second, it enables the effects of pitch-flat/edge coupling to 
be systematically introduced. 

Rigid- Body Feathering Degree-of- Freedom 


For this purpose the total pitch angle now includes the root feathering 
angle, 9^, and is expressible as; 


9 
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( 62 a) 





where 


% 

''ntm+i 



(62b) 


Qii = 9o (root torsion deflection) (62c) 

t'NTM+l ” 


A separate rigid-body torsion equation is written by a generalized appli- 
cation of the Galerkin technique wherein the torsion equation is multiplied by 
the pseudo-mode and integrated; this effectively calculates the torsion moment 
resisted by the equivalent elastic root torsion spring, K02: 


/ [distributed torsion loodingjdr = 

(GJ)pa 


»FB 


1 ■ "J 


(63a) 


where : 


= torsional stiffness of flex-beam = ( 63 b) 

K^s equivalent root torsion spring due to control system flexibility, etc. 


In addition to incorporating Eq. (63a) into the total set of dynamic 
equations, addltioneG. terms must be added to the equations for the J normal 
torsion modes : 
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T^Jdistributed torsion loading] dr = J ^.Qax 




(Gj) 


FB 


9 . 


’FB ®FB 


(64) 


Furthermore, appropriate 9 r (and time derivative) terms must be added to the 
load distributions (wherever ® dependent terms appear). 


Inclusion of Pitch-Flat/Edge Coupling 


Automatic pitch change coupling effects can also be conveniently included 
in the equations of motion by \ise of the pseudo-torsion mode. First, however, 
the kinematics of the push-rod to blade attachment point must be considered. 
Let Az be the upward "5" coordinate system displacement of this attachment 
point due to blade deflections with the push-rod momentarily disconnected. 
These deflections are shown in Fig. 5 for a blade section at the attachment 
point station: 

. Z5 



Figure 5» - Geometric Jteatxires of Automatic Pitch Change 
■ due to Blade Deflection. 
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When the push-rod is reattach'ed, so that the attachment point is restored 
to its original position, the blade will have rotated through an incremental 
pitch angle (-A0 as shown in the figure). Using simple trigonometry the fol- 
lowing equation governing A0 and Az can be written: 


sin^Opp- sin(0OpR+^^) 

which upon making the small angle 


= 

VlOpR 

assumption on A0 becomes: 


(65) 


A0 



Az 

VlOpR 


( 66 ) 


It should be noted that the Az deflections can result, depending on rotor con- 
figuration, from flatwise bending, edgewise bending, flapping and/or lagging. 
Thus, A0 can generally be written as: 


A 0 = 0 Wjqwi + \qv^+%)3+ei5S ( 67 ) 

and the deflection of the reference point of the blade attachment point section 
(point 0, in Fig. 5 ) can similarly be written as: 


( 68 ) 



•Hie above development can then be incorjKsrated into a comprehensive 
mocteling of pitch- flat/edge coupling effects by building iipon the following 
two considerations. First, the incremental pitch angle, d0, can be incorporated 
Into the total pitch angle and twist rate, again using the pseudo-torsion mode: 

NTM 

0 = E 

js| J J 


NTM 

»■ = E »o + 8^ (70) 

j S I J J 


Hence, the incremental pitch angle generates aerocSynamic and inertia loads 
pjropprtional to those which would be generated by the root deflection angle, 
0j^, whose dynamics are described by Eq, (63a), 

Ttie second consideration is that in the process of equilibrating the 
blade in torsion and maintaining the push-rod attachment point at its level 
position (as indicated above in Fig, 5) a push- rod force is generated equal 
in magnitude to the applied torsion moments on the blade divided by the push- 
rod offset (Slstance, but directed in the negative direction: 





(71) 


where the torsion moment is obtained from Eq; ( 63 a) and given by: 


M 


[distributed torsion loading] dr 


- +A0 + r.(rj) +«3(rj)J + -s^e^ 


'FB 


(72) 
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Ibe push- rod forced .> is then a (negative) concentrated force acting on 
the virtual displacement iMs displacement Is derived from Eg,. (68) and la 
escpresslble by the following equation 

8zg = AZjSqy^. + AZ|j8qy|^ + Az^8/9 + AZg8(6) (73) 


Bierefore, combining Eqs. (71) and (72) -and using Eg. ( 67 ) j the bending «gua> 
tlons and those for flapping and lagging are modified as follows: 

l*th flatwise bending. Eg. (24) ; 




AW: )( inertia loads) + . . . } dr + M- 

' ■' i 5pR 


f {(Xyy.- AW. )( aerodynamic loads) + ...}dr 


. s O 


(7*») j 

I 


k*th edgewise bending. Eg. (2g); 


J inertia loads) + • . • } dr + 

+ AV|^)( aerodynamic loads) + • • }dr + . . . * O 


rigid flapping. Eg. (31) ; 




( 76 ) 
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rigid lead-lagging. Eg., (32) ; 



+ q2^)dr + 6^M 


*5 PR 


+ . 


= 0 


(77) 


The inertia matrix, Ca], developed in an above subsection can then be 
modified, as a result of Eq,s. (69) and (73) thru (77 )> to the following general 
forms; 


Apo = ApQ+ 9 pA^^Q + (78a) 


Ap^. = Ap^ + 9 pAgi ( 78 b) 

J J " J 


A^.o = ^djO + ®oA0j0^ ( 78 c) 

where P and Q denote any of the subscripts, Wj^, Vj^, p or 5, and i^ere Aj^j denot: 
any of the basic inertia matrix components given in Eq,s. (46a) thrru (4-6n;. !IM- 
technique for the introduction of a finite push- rod force in the bending, flap- 
ping and lagging equations is appropriate only for the case wherein the root 
torsion retention stiffness is asstimed infinite, no rigid feathering degree-of- 
freedom exists and the incremental pitch angle is automatic. For cases where 
a finite root stiffness exists the push-rod force is expressible more conve- 
niently using deflection dependent spring forces . 

The above somewhat abbreviated development is sufficient to write all the 
explicit terms involving automatic pitch change due to blade bending, flapping, 
and/or lead- lagging. For the bearingless rotor with a cantilevered torque tube 
however, this development must be expanded to include the finite effects of 
flexbeam bending. This bending gives rise to the low stiffness torsion (wob- 
ble) mode whose equations are developed in the next subsection. 


AEROEIASTIC SIMULATION OF BLADE TORSION WOBBLE MODE 


The cantilevered torgue tube configuration bearingless rotor cong>rises a 
relatively siniple mechanical system. The torque tube is attached at its in- 
board end to the push-rod and at its outboard end to the flexbeam and outer 
blade portions of the blade by means of a cantilever mount (see Pig. 6), 
Typically, the skin of the outer blade portion would continue inboard from the 
flnbeam-^ torque tube Juncture to form the hollow torque tube. As shown in 
Pig, 6, this configuration is characterized by a relatively long shear load 
path for the push-rod loads. 

Three observations can be drawn ftom this figrore. First, because of the 
combination of flexbeam flexibility in both torsion and bending, the blade 
possesses a rigid body feathering degree-of-freedom even with infinite control 
system in5)edance. This torsion degree-of-freedom inherently involves signifi- 
cant vertical or flatwise motion and is referred, herein, as the ”wobble” mode. 
Second, because flexbeam bending plays a primary role in defining the restor- 
ing moment for this mode, an analytic calculation of this flexibility is re- 
quired. Third, as a result of this flexbeam bending flexibility the conven- 
tional automatic pitch change coupling described in the previous subsection is 
no longer completely applicable. In the present case, the effect of bending 
modal deflection is to develop a push-rod load which acts as an applied torque 
to the wobble mode, which in turn has finite impedance characteristics. Thus, 
for the cantilevered torque tube configiu:*ation effective pitch- flat/edge 
couplings may, in general, be dynamically amplified and phase lagged. The 
development of the equations governing the wobble mode presented herein draws 
upon the pseudo- torsion mode formulations of the previo\is subsection. 
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Figure $• - Flctorl6LL Representation of Wobble Mode and Load Bath for Shears due to Bush-Rod Loa4> 



Flatwise Flexibility Characteristics of 
Flexbeam - Ceuitilevered Torque Tube Assembly 

Ihe reqTiired quantitative description of the flexbeam bending flexibility 
is the flatwise deflection of the inboard end of the cantilever mounted torque- 
tube per unit load applied at the same point. Ihe assumptions made for this 
calculation are 1) that the torque tube is rigid in bending, 2) that the 
effects of flexbeam twist angle can be neglected, 3) that over the flexbeam • 
span the tension is constant and equal to that due to totHl blade mass (l.e., 
the flexbeam is massless) and 4) that the junction point is at the midpoint 
of a blade segment. The elastic problem so fomnilated is solved by means of 
transfer matrix techniques using the distributed (lumped) mass smd flat^se 
stiffness properties assmed for calculating the normal flatwise ben^ng modes. 
Figure 7 shows the features of the mathematical model of the total blade. It 
is to be noted, furthermore, that all the blade segments are in a centrifugal 
force field. 


P-1 lb. 



Figure 7. - Features of Mathematical Msdel for Calculating Flatwise Flexibility 
Characteristics of Inboaz^ End of Cantilevered Torque Tube. 
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Since the junction point is assumed to act at the midpoint of a segment 
and since the centrifugal force is a strong (hut known) function of span, each 
blade segment is further broken down into two semi-segments, each of which is 
then assiamed to have a constant (average) tension. Ihen, for each semi-segment, 
the inboard loads and deflections, expressible as a state vector, can be related 
to the outboard loads and deflections by means of an appropriately calculated 
transfer matrix; 


(s 

1 1 

fs) 

) M 1 

. - TtI 

M ( 



z ( 

(z'J 

1 1 

[r) 


k = l, 2,..., 2(NSEG) 


(79) 


where : 


N 


L+ -l)c| [i +(C 0Sh;8L-l)C 

+ (C0Sh/SL-l)C2 


Tf.sinh^c, 


+ Cg/8sinh/3L 


(80a) 


-CgCOsh/SL+Cj 


±(cosh^L-l) I 


--“(cosh)8L-l) 


- -^sinh)8L 


sinh^SL 

“7“ 

cosh/3L 
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Qjie various terms are evaluated at the k*th semi -segment and defined as foUovs: 


T = average tension = i/2(Tj+Tq) 


(80b) 


= yx/EI 


L = I/2AX 


C, = =- To +mL‘ 


^2 L 20/] 


* Cl 


nr» / _L . I \ 
■ T \ 12 L 30/ 

- I ' \ 

T V '2 L "^20/ 


and where: 



k = l corresponds to the 
blade tip semi-segment 


( 80 c) 

( 80 d) 

(80e) 

( 80 f) 

( 80 g) 

( 80 h) 


Cascade matrix multiplication, together with the Introduction of the 
additional shear and moment at the Junction point (arising from the unit load, 
P) and the imposition of boundary conditions (geometric at the root and loading 
at the tip) leads to a matrix equation of the following form: 
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iHiis matrix equation constitutes a solvable set of four equations in four 
urtoowns. A pstrtial cascade multiplication of transfer matrices, T^j relat- 
ing the jTinction state vector to the tip state vector, is then used to ob- 
tain the deflection and slope at the junction point; 


( 82 ) 


Ihe final required flexibility is obtainable from these two quantities as 
follows ; 



Z = 


Zj-SttZj + 



( 83 ) 


where : 


Sipj =» length of torque tube 



flatwise angular spring restraining the torque- 
tube to the blade spar ( flexbeam ) 


Elastic Rigid Feathering Torsional Restraint 


The combined effects of flexbeam bending, flexbeam torsion and control 
system flexibility can be determined by examining and summing the contributions 
to the vertical deflection of the push- rod attachment point. As shown in Fig. 
8, the attachment point is allowed to deflect sequentially due to the effects 
of blade bending (Az), of flexbeam bending due to the application of the push- 
rod load (zg), of flexbeam torsion (6 r), and finally of control- system flex- 
ibility (Z3). 

The final attachment point deflection is thus given by the following 
expressions : 




(84) 




Figure 8. - Contributions to Vertic€LL Deflection of Push-Rod Attachment Point, 
Cantilevered Torque Tube. 


Oils deflection together -with the actueJ. control system stiffness. Kg*, 
defines the push-rod load; ^ 

K K 

fpi. = = -jF^ (85) 

^•OpR ’^'OpR 

idler e: 

K0^ ■ actual control system stiffness. 

Recognizing that the Z2 deflection is Itself proportional to the push-rod 

load: 
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( 86 ) 


Z2 = Z*Fpp 

enables the following eaqpressions for the push- rod force and root moment to be 
written: 


Ka /y2 

F = 3 ^lOpR 

*^3 'lOpR 


(^z + ViOpr^r) = 




^2 
lOpR 


(^2+9,0pr^r) 


(87a) 


(87b) 


where the effective root spring, introduced in the previotis subsection, KQ^f is, 
now given by: 




= 


l + K. 


,2/y.o 


PR 


( 88 ) 


lynamic Equation for Wobble Mode 


As was observed above in this subsection the wobble mode is essenti ally a 
rigid body feathering mode; hence, the developments of the previous subsection 
can be drawn upon. First, the quantity A 0 appearing in Eq. (87b) for this 
cantilevered torque tube application, must now be considered a normalization 
of the flatwise displacement of the attachment point rather than, as in the 
previous subsection, an automatic pitch change. Thus, the total pitch angle 
and twist do not contain the quantity A 9 and Eq. (6la) and (6lb) become the 
appropriate expressions for the total pitch angle and twist, respectively. 
Using the root restraint moment, as defined by Eq. (l. 87 b) the dynamic 
equation for the wobble mode then becomes a modification of the rigid feather- 
ing Eq. (1.63a) given in the previous subsection: 


/ )^^B[‘^'Stributed torsion loodingjdr = 


FB °fb 




(89) 
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Coniparison of this equation with Eq. (63a) shows a coupling of the wobble 
mode with blade bending through th^ A 6 term. It is reasonable to e^^ect a 
similar coupling of the bending equations with rigid body torsion through a 
similEur tenn. Indeed, the push-rod force, given by Eq. (87a) can be directly 
incorporated in the flatwise and edgewise bending equations to give the ex- 
pected symmetrical elastic coupling: 


/ {t^w.~^Wj)(inertia loads) + ...}dr + 

+ J |(y^,-AW.)(oerodynamic loads) + ...}dr+... = 0 


(90) 



AV,j)(inertia loads) + . . . } dr + - a9) 


+ y*|(yy^-AV^)( aerodynamic loads) + . . .}dr + . . 


= 0 


(91) 


The development of this subsection has been directed to a type of beaaring- 
less rotor configuration which is characterized by a torque tube which is 
relatively stiff in bending. This high degree of torque tube bending stlffhess 
is generally reqtilred because of the inherently long load path for the push-rod 
shear loads. In the next subsection an alternate type of torque tube is con- 
sidered wherein the torque tube bending stiffness is not generally high. 3 Ms 
characteristic leads to a degree of bending redundancy in the flexbeam (torque 
tube elastic system) and an appropriate analysis technique is required. 
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EEDUNDMT ANALYSIS OF FLEXBEAM - OJORqUE 
TUBE ASSEMBLIES WITH SNUBBER RESTRAINTS 


Basically, bearingless rotor systems are characterized by an inboard span- 
wise section comprised of two highly specialized and contrasting structural 
elements (see Fig. 9)* first is the innermost portion of the blade spar 

and is referred to as the flexbeam . Its function is to provide the inboard 
bending stiffness appropriate to hingeless rotor blades, while at the same time 
being very soft in torsion to allow the outer portion of the blade to be feathv- 
ered. The second structural element is typically a cylindrically formed shell 
relatively stiff in torsion which encloses the flexbeam and is referred to as 
the torque tube. The structviral functions of the torque tube are to provide the 
blade torsion system with sufficient general torsion stiffness inboard of the 
junction and to transmit those torques to the flexbeam which are needed to 
produce the input control angles. The inboard end of the flexbeam is attached 
directly to the hub whereas the inboard end of the torque tube attaches to the 
control push-rod eind to the flexbeam . More precisely, the inboard end of the 
torque tube is restrained in torsion by the push-rod and control system, and 
in shear by means of an (effective) pinned point (elastomeric snvibber or equiv-. 
alent) attached to the flexbeam . The cantilevered torque tube described in an 
earlier -section is restrained- both in torque and shear by means of the push-rod. 
The red\indant analysis of this latter configuration requires a much more exten- 
sive development than is presented herein for pinned or snubbed configurations. 
Because of the aarrent viability of the snubbed configuration, and to a lesser 
extent, the additional complications of analyzing the cantilevered configura- 
tion, the present development is restricted to only pinned or snubbed configura- 
tions. The two structural elements are attached (both to each other and to the 
remaining outer portion of the blade) at their common outboard end point, 
referred to alternatively as the junction or juncture. 

For typical bearingless rotor configurations the torque tube and outer 
blade portion skin would be integrally formed (for practical fabrication reasons). 
Hence, while the torque tube would nominally experience zero bending moments at 
its inboard end (due to the effective pinned joint), it would necessarily sup- 
port bending loads at its outboard end. Eqviilibrating the internal loads at 
the junction point arising from the load transmissibilities of the flexbeam , 
torque tube and outer blade portions of the blade defines a doubly red\mdant 
analysis. Over the flexbeam torque tube portion of the blade, dual load paths 
are defined in both bending and torsion, and an appropriate solution must con- 
sequently require that the elastic deflections of the jiinction point, both in 
bending and torsion of the three structural elements, are consistent. The 
remaining portions of this section develop the detailed mathematical formula- 
tion and solution to this redundant analysis problem. For conceptual clarifi- 
cation the primary bending flexibility system is considered to be the flexbeam 
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whereas the primary torsion system is considered to be the torque tube. The 
subsections to follow include; first, a development of the bending stiff- 
ness characteristics of the torque tube; second, a development of the 
torsional stiffness characteristics of the flexbeam arising from nonlinear 
loadings; third, a method for estimating the internal loads immediately out- 
board of the Junction; and finally, a mathematical synthesis of these elements 
to effect the complete solution to the redundant smalysis. 


Bending Stiffness Characteristics of the Torque Tube 

The appropriate elastic bending description of the torque tube is a 
stiffness matrix expressing the inboard shears, the outboard shears and the 
outboard moments as linear combinations of the relative deflections of the 
ends of the torque tiobe. To this end the following specific assumptions are 
made: 


1. The torque tube is a beam whose section properties have generally 
nonuniform spanwise variation. 

2. The torque tube has negligible twist but is inclined by the pitch 
angle of the Junction point. Hence, flatwise and edgewise stiffnesses 
generally couple the in-plane and out-of -plane bending characteristics. 

3. The spanwise in-plane and out-of -plane load distributions over the 
torque-t\ibe span act directly upon the flexbeam ; i.e., the torque tube 
is transparent to these loadings. 

4. The radial loadings due to torque tube mass (centrifugal forces) con- 
tribute to the stiffening of the torque tube in bending. 

5. In addition to having distributed bending stiffnesses, the torque 
tube is elastically restrained at its outboard end with finite angular 
springs (defined in torque tube flatwise and edgewise orientations) and 

at the inboard end with finite lineal shear springs (defined alternatively 
either in torque tube flatwise and edgewise, or in hub out -of -plane and 
in-plane orientations ) . 

6. The torque tube is supported radially either at its inboard end 
(torque tube in' tension) or at its outboard end (torque tube in compression). 

Assuming the torque tube to be a beam allows a transfer matrix solution 
similar to that employed in the previous section for calciilating the flatwise 
flexibility of the flexbeam . Analysis of a beam segment with constsuit pro- 
perties forms the basis for evaluating each component transfer matrix and 
follows from a consideration of such a beam segment shown below in Fig. 10. 
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Figure 10. - Loads and Deflections of a Beam Element 
with Constant Section Properties. 


Use of the standard beam theory requires the bending equilibrium equation 
for the k'th element to be written as; 




(92) 


where the average tension is: 


\ Xok> (93> 

As shown in Fig. 10, Tj^ represents the average tension over the beam 
element. Note that for torque tubes in compression would be a negative 
number. The bovindary conditions appropriate to E<^. (92) are: 


63 


L 



( 94 ) 
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* • *0,,. *'• *0,, 0lX.L|,| 


For torque tube elements in tension the above formxilation duplicates that 
given in the previous section and the resulting transfer matrix, as given by 
Eqs. (80a,b,c, and d) is applicable to this case. For the case of a torque 
tube in compression, the quantity, 3, is then defined as: 


^ » y-T/EI » -/C/EI 


(95) 


and the resulting component transfer matrix is given as: 

n I 0 0 0 



[i +(i -cosiSUdi 
+ /3sin/3Ld2] 

,-r sin/9L ^ 

0 

-dgC0S/9L + ds] 

(96a) 

1 /, sin/3L\ 

C r } 

id- cosy3D 
C 

I sln/3L 

0 


^ (1 - cos/SD 
C 

sin/3L 

0 cos/3l 

k 

d, = y |-Cq + 



(96b) 

d - / Xi 



(96c) 

*^3 ' — Vl2L 

io) 


(96d) 
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As in the previous section, the total (uncoupled) stiffness matrices, 
for both flatwise and edgewise bendli^ . are obtained by a cascade multipli- 
cation of the component transfer matrices as Indicated below: 

[r].n[Tj (97) 


The uncoupled e<3gewlse and flatwise stiffness matrices are defined In the 
"6” coordinate system and are denoted as: 





(98a) 


(98b) 


where the subscripts (S) and (j) denote, respectively, the snubber and Juncture 
ends of the torque tube. Eqs. (98a and b) must then be coupled by the effects 
of Juncture mcsment springs, snubber shecur springs and torque tube pitch angle, 
6j. The Juncture moment retention springs are defined by the following simple 
relationships : 


(Mj 


6^J 



(99a) 
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F 


(99b) 
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similarly, for Inboard stiubber shear springs aligned In the torque ttibe 
flatwise and edgewise directions ("6" coordinate system), the following rela- 
tionships hold: 





^6- ^ 


'FB 


(100a) 



- Z* 


*TT 


)*K, 


•FB 


(100b) 


Noting that and are both zero (snxibber transmits shear only), 

Eqs. (99a 8ind b and 100a and b) can then be combined with the definition for 
the uncoupled stiffness matrices ( 98 a and b) to give the following matrix 
equations . 
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In both of these matrix equations the respective premultlpllcatlve 
matrices on the left hand sides can be inverted to give modified stiffness 
matrices denoted in the following abbreviated form: 


1 







(102a) 
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(102b) 




These equations can be represented even more compactly using the follow- 
ing notation: 



^ 1 ^ 

O |Rp 


AZe 


(103) 


Beca\ise both edgewise and flatwise transfer matrices were derived 
"uncoupled" (i.e., each assumes zero deflection for the other), the loads 
given by Eq.. (103) are properly interpreted to lie in the planes defined by 
the torque tube sections at the snubber and at the juncture. These pletnes are 
oriented by rotations about the y^ and z^ (or yg and zg) axes by amounts deter- 
mined by the blade slopes at the snubber and juncture points. For what 
follows, the subscript "7" is given to the moments so defined by the torque 
tube stiffness matrix with the understanding that they lie in the plane defined 
by the torque tube section at the jvincture. 

Transformation to the "5" coordinate system using trigonometric resolu- 
tions (rotations through the angle, 6j) on the elastic load vectors, Ly^ auid 
Lgg, and the deflection vectors and AZg then yield the following coupled 

form: 

i LygJ |“R£C0S^ej + RpSin^8j I (RE-R p)sinejCOsej ~] (ay^) 

Lrg ' [(Rp-Rp)sinejcosej | RESin^8j + RpC0S^ejJ ( ,(104) 

or, in abbreviated notation: 

M” ■ MN,. 

where the (d) superscript denotes loe^ds due to the (A^)^ def4.ectlons 
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Should the inboard snubber springs be aligned in the hub axis system 
rather than in the torque tube flatwise, edgewise system, then Kgg“l and 
in Eqs, (lOla and b) are zero and the following description of the 
snubber elasticity, must be used: 


<V*' ‘'’5— 


•TT '»>Fe 


(S, J_ = (2g -2. )-K- 

*5S SspB 


(I06a) 

(I06b) 


These linear relationships are then inserted into the elastic description as 
follows : 

W) 


1 

1 

JP 

1 ffi 1. 

0 

-S|4‘<l25 



\ 

I 3 

\ 

"S24Ks^ 

*■5 

-534^5 

-S44»<1z5 

0 

y5 





-S54Ki,^ 

0 

0 

-S64k1j2jj 

.1 

\ 



-S74Ks,^ 

-S 84 KI 2 JJ 

I 3 

\J 


\ 


-ys 


> ■ 


\ J. 


(107) 


TT 


where Kg and Kg^ are the snubber springs in the in-plane and out-of-plaine 

direction!, respectively. Eq. (10?) can then be inverted to yield the desired 
form, similar to Eq. (105): 
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(108) 


where the subscript (FB) on the right-hand side indicates deflections of the 
flexbeam at points where the torque tube and/or torque tube retention springs 
are attached. Note that the torque tube inboeurd slopes are Included in the 
load vector first to indicate that they, like the loads, result from the spe- 
cified deflections on the RHS of this eq'uatlon, and second, to show where in 
the analysis these quantities are formed when needed for pitch-flat/edge 
coTjpling calculations. In practice, these quantities, once formed €U 7 e parti- 
tioned out and handled separately. In subsequent development they will be 
omitted from the load vector, for clarity. 

The total internal loads acting at the ends of the torque tube consist 
of those elastically produced by the deflections, given by the above equation, 
and those loads arising from blade root torques and the resulting pvish-rod 
loads. The torque tube and loads due to push-rod loads are already implicitly 
Included in the flatwise and edgewise bending excitation equations by reason 
of the particular method for including pltch-flat/edge coiipllng discussed in 
an above section. These loads, however, must be explicitly evaluated to 
Implement the torsional redundancy portion of the redvindant analysis. The 
z«mainlng portion of this subsection discusses those derivations needed to 
calculate the pltch-flat/edge coupling when a flexible torque tube is en^loyed, 
how £q. (108) for elastic flexibility is used to account for the bending redun- 
dancy portion of the redundant analysis, and how the torque tube loads due to 
push-rod force can be estimated. 

As was described in an above section and shown in Fig. 6 , pitch-flat/ 
edge coupling effects are calctilated from the geometry of the pitch-horn and 
from knowledge of the vinconstralned vertical deflection of the attachment 
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point, Lz, for unit deflections of the various modes. With the above einalyt- 
ical description of torque tube bending flexibility this can be readily accom- 
plished. The Az -deflections per modal deflection for this case can be obtained 
from the last row of Eq. (108) for and from (z5g)j^ which is either 

eqvial to or, if the snubber spring rates are finite, calculable from 

the linear relationships given by Eqs. (lOOb) through (105) and (l06b) through 
(108) using routine algebraic manipulations. After these manipulations are 
performed, (z5g);piji and ( and, by linear combination, Az, can be ex- 

pressed as a linear combination of the flexbeam deflections due to unit modal 
deflections ; 


A2.(Z5^)^-Mrp„-r5)(^5^V^. [qttJ 




(109) 


Zgj /per unit modal deflection 


The above equation for Az is then used with the derivation in Section V to 
obtain the values of modal pitch-flat/edge coupling. 


The bending redundancy portion of the redundant analysis is acconplished 
by treating the snubber and Juncture loads, defined by Eq. (IO8) as concen- 
trated loads for direct inclusion in Eqs. (24) and (25) as is provided for in 
these response equations. Due to the linear character of Eq. (IO8), the 
bending redundancy can be included in the eigensolution as well as the time- 
history solution. 


Calculation of the external torque tube and loads resisting the push-rod 
force is accomplished by equivalencing the virtual work due to the push-rod 
load with that due to the equilibrating torque tube and loads. To this end 
Eqs. (109) can be used in the following manner; 

Sv = fpR'SAz = Fpp -S (110) 
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Therefore, since -the con^onent virtual displacements are Independent, 
the externeil end loads needed to equilibrate the push-rod force, Fpjj, are 
given by the following: 


^TT. 


-q 


(PR) 


*• •'TT 


TTg 


PR 


( 111 ) 


which, when combined with Eq. (108), yields the total internal loads at the 
ends of the torque tube: 


r 1 r f V™’ 

{'■*}„ ' {‘•'Irr {'-'It 


'TT 


'TT 


= ['.] 


>'»s 

y«» . 
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‘^TT, 

^TTe 

^TT3 

I -^TT^ 

'’tT,: 


-q 


TTe 


PR 


( 112 ) 


In addition to the loads given by the above equation, the radial shear 
at the Juncture, is given by the following: 


(S 




torque tube in tension 



torque tube in compression 


(113) 


Note that the axial torsion moment at the outboard end,. (Mxy )^^, includes 
the effect of the flexbeam resisting moment and, hence, mus^ be solved for 
using the redundant analysis. The following siibsection formulates the analysis 
stiffness matrices required to define the detailed torsional stiffness of the 
flexbeam . 
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Torsional Stiffness Characteristics of the Flexbeam 

Calculation of an appropriate stiffness matrix statement for the flexi- 
bility of the flexbeam , requires that two special characteristics of the flex- 
beam be accounted for; First, the coni)osite material flexbeam is a plate- 
like beam whose cross section has a typical aspect ratio of order of magnitude 
10; furthermore although the flexbeam is relatively soft in torsion (low GJ), 
it retains a relatively high bending modulus (high El). Hence, it is reason- 
able to expect that the appropriate differential equation governing the 
torsional elastic characteristics should contain the torsional stiffening 
effects dvie to plate bending. Second, the B]_ and B2 constants were retained 
in the general development of the blade torsion equation (Eq. 30) principally 
for applicability to non-helicopter rotor applications (i.e. propellers and 
wind turbines). However, the effect of these constants on the total stiff- 
nesses of composite bearingless (helicopter) rotor flexbeams is considered 
to be. negligible . Indeed, assuming a large section aspect ratio (^ ^ 10), 
the ratio of twist squared stiffening to St. Venant (GJ) stiffening can be 
approximated by the following torsion stiffness ratio: 


For comparable section and planform aspect ratios (c/t and l/c, respec- 
tively) and even for an (e/g) ratio of 25, typical of unidirectional cstrbon- 
epoxy sections, the above torsion stiffness ratio reduces to approximately 
O.UAG^. The total flexbeam twist angle, A9, then must be in excess of 25 
degrees to achieve a torsion stiffness ratio of 0.1, which is considered 
justifiably small for present purposes. Third, because of the length of the 
flexbeam and the nature of the bending loads and deflections at its outboard 
end (resulting from torque tube and outer blade generated internal loads), the 
torsion moment over the flexbeam span is generally variable and contains 
components due to all six concentrated loads at the juncture. These three con- 
siderations thereby define the fo3J.owing differential equation for a torsional 
element with constant section properties: 


DC^ 

I2R^ 


9 + 


-^(Gj+TkJ)^' = 


*9 


(il4) 
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'Wherein boundary conditions on 9 euid 9* must be specified at the inner and 
outer ends of each torsion element. At the ends of the total length flex- 
beam (inboard end (k = 1), and the junctTare end (k = K)), 'the appropriate 
boundary conditions are: 


0(0) = 0'(O) = ©Vj) = 0 

0(rj) » ®j"®Bp0= 


(115) 


Before solving Eq. (ll4) let us first examine the structure of 'the right- 
hand side. Figure 11 below shows the origins of the spanwise variable Internal 
torque which accrue from combinations of loading and deflection at the out- 
board (;Juncture) end of the fl exbeam . 


*^2 



UOADS APPLIED TO 
FLEXBEAM BY' 
TORQUE TUBE AND 
OUTER PORTION 
OF THE BLADE 


Figure 11* Pictorial View of Flexbeam Showing Deflections cmd Applied Loads 
Producing Spanwise Variable Torsion Moment. 
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Eq, (27) can be readily adapted to’ the present case irtierein the loads are con- 
centrated at a point; the following e^^ression for the spanwise variable 
torsion moment, c&o then be written as: 


M, (r) 


= 25)75 - (Y5j- 75)25] Sxg- [z5j"Z5 " <S- D25] Sy^ 

+ [v^5 - (S -r)7s] + 75M75+ z's %+ % 


(116) 


since the span of the flexbeam Is short relative to the total blade spaui, 
the In-plane and out-of -plane deflections, anda^, respectively, can be 
approximated by polynomials defined by the deflections and slopes at the Junc- 
ture. With the definition of a flexbeam spanwise variable: 


1 • r/S. 


(117) 


the deflection (and corresponding slope) in either of the two directions can 
be expressed by the fo 3 J.owing general expressions: 


(75 or 25) = f = (3fj-§f'j)i7^r(2fj-Sfj)-»7^ 


(118a) 


(75 or 25) = f'= i[(6fj-2Sfj)i7“(6fj-3gfj)77^] 


(U8b) 


where : 


S = S/R 


(118c) 
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By using Eg.s, (U8a and b). Eg,. (U6) can be expressed as a polynomial 
in the flexbeam spamrise variable, 11; 


MxgW) = + +(l-6,7+V-477^(Sz5y5j- 

+ 2S(^-2'7^+i7’)(Szj.y5j-SyjjZ5j) + -|(^-’7^)(My5y5j+ 


Thus, the right-hand side of Eg. (U4) is seen to be a polynomial of 
fourth degree in the spanwise variable 7). Since the flexbeam , like the torgue 
tube, will have section properties variable with span, a transfer matrix solu- 
tion using span segments must be employed. For any constant section spanwise 
segment the appropriate general solution to Eg. (Il4) becomes: 

e = ce^’^ + De^’' + i HjV (120) 

i =0 


where C, D and are constants to be determined from the boundary conditions 
eind from the polynomial coefficients of Eg. (119), and where the section 
characteristic constant, X, is given by: 



(6J, +Tk|) 


( 121 ) 


Since all the botindary conditions given by (115) csui be used only after 
all the component transfer matrices are cascaded, one of the constants, 

Hq specifically, must be treated as a state variable (like 9 suid 0’) and Varied 
by each component transfer matrix. Eg. (120) can then be used to form an 
intermediate (partial) transfer matrix eguation: 
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d 

/ 

e 

Ho 


coshXi --j^sinhXi 
-XsinhXi coshXi 

0 0 


(i-coshXjt) 

XsinhXii 

I 



/ 




III \ M sinhXJt 

+i’7o X — 

'■^0 ^OShXA 

0 


( 122 ) 


•where, far each k'th segment, •the constants = 1, 2, ..., 5) are com- 

pletely de-termined by the section properties and the polynomial coefficients 
of Eq., ( 119 ). The length of the k'th spanwise segment, is conveniently 
taken to be the blade break-up segment length, Axj^, divided by s. More speci- 
fically, using the following notations: 








B„. ^(GJ+Tk®)|, 


( 123 ) 


the explicit expressions for can be written as: 


Bk 

<|m 

CVI 

+ 

<l|CD 

CM 

+ 

hP 

(124a) 

11 

CSJ 

X 


(124b) 


>Tk+4(|.)^X,] 

(124c) 

H . -i- 

H4k 48 

k ’ ” 5 I> " 5 ^ '''4 

(I24d & e) 


Equation (122) is deemed a partial transfer matrix because it does not 
yet include the transfer of the Hq constant from one segment to the next; 
one additional relationship is thvis required to comple-te the transfer matrix 
formulation. For this purpose, Eq. (ll4) can be Integrated across an in'ter- 
segment junction point, (T) = 'Pjj. ' 
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[-A9"+ - [-A/+ . 0 (125) 


k,k+i 'k,k+i 


Since the torsion moment is a smooth continuous function over the entire 
span the integral in the above equation must equal zero. Therefore the quantity 
(-A0'*+B0) must be continuous across the inter-segment junction point. This 
relationship results in the following additional required equation: 

[-A9%B9] ■ f:Hik[B^4^^,,-A^i(i-l)^«^,]: , 


®k+i ^O, k+i 




^k+i’^k.k+i ^k+i'^' '^ik, k+i] 


(126) 


Combining Eqs. (122) and (126) yields the following final expression for the 
complete transfer matrix equation: 

0 

\ cosh(XA)^ -- jJj-sinh(XA)^ [i- cosh(XJt) J ^ 

\ ® / " -X,^sinh(XJl),^ cosh(XX),^ | ^ 1 


®k+i ^'®k 


•»7j - ■»7 q cosh XI + -^ Vq' sinh xi 


( I ) + 17^^ Xsinh Xi - ( i ' cosh Xi 


I, {■<* (^ V 


(127) 


I - coshXi 
XsinhXi 
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Note that the constants, (i = 1, 2,.., .,5) in the above equation are com- 
pletely determined by Eqs. (124a throxigh e). This transfer matrix equation 
for the k’th segment can then be s\iitably cascaded to yield a final matrix 
equation, relating the state variables at the flexbeam root to those at 
the Juncture: 



The E and F matrices appearing in this equation are the result of appropriate 
cascading and, hence, are not sufficiently siniple to express explicitly herein. 
They are calculable, however, in a straightforward manner using routine matrix 
algebra. 

Application of the boundary conditions, Eq. (115) together with the 
elimination of and (H^)j^ from the above equation set results in the 

following expression for the elastic torsion deflection of the flexbeam at 
the Juncture: 


1®||®23“®2I®I3) i^o 


’’’i^^2,i + l®l3 ^l,i + l®23^ 


(129^) 


where ej^j and f.j^ are the elements of the E and F matrices, respectively, and 
where from Eq. (119) : 
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9 
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1 

0 

0 

0 

0 
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SzgYsj-SygZsj 

SzgYsj- SygZsj 




^5^5, 


+ M, Zk 
I ‘5 ®J 


(129b) 
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Equations (l29a and b) can be cdhbined using the following compact notation: 


A6j = %x^*5^y5jZ'5j~ ZSjVSjVb 

+ Syg Z 5 j)pg + ©s/SzgYgj - Syg Z 5 j)pg (l 30 a) 

+ »Mo^'^y5y5j+ >^Z5Z5 o)fB + 


or alternatively: 


A0J 


= ^SxSx5pg+ ®Mx>^X5pg+ S^ySpB 


( 130 b) 


where the various stiffness coefficients in Eq. (l30b) are deflection dependent. 

Thus, the elastic torsion deflection, A0j, is seen to result not only from 
the usual radially oriented feathering torque, Mx^, but from various nonlinear 
combinations of loads and deflection. The following subsection presents those 
calculations needed to estimate the internal loads just outboard of the junc- 
ture which must be equilibrated with the torque tube and flexbeam loads at 
the juncture. 


Estimation of Internal Blade Shears and Moments 
in Blade Adjacent to the Juncture 

Two basic me"thods exist for estimating the instantaneous internal loads 
in the outer portion of the blade: mode deflection and force-integration. A 

comprehensive discussion of these two methods and their accuracy characteris- 
tics is contained in Ref. 8 and will be omitted herein. Briefly, the former 
is relatively simple to implement but requires smooth speinwise load distribu- 
tions for accurate estimates with a small number of modal variables; the 
latter is, prima facie, a more complicated method to implement but estimates 
the loads with superior accuracy. The latter method (force integration) was 
chosen because, for most practical helicopter applications, it is significantly 
more accurate, and because the disadvantage by reason of complexity is molli- 
fied by certain available features implicit in the general analysis solution. 
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Because the torsional redundance portion of the redundant analysis Itself 
entails significant nonllnearltles, th6 redundant analysis effects are most 
meaningfully Included oiily In the total nonlinear formulation herein the 
equations are solved by numerical Integration to obtain time -histories . Exam- 
ination of the various differential equations formulated for the aeroelastlc 
responses, Eqs, (24), (25), (30), (31), and (32) reveals that the Inertia and 
aerodynamic load distributions are used explicitly. That Is, with the excep- 
tion of the double (hondlmenslonal) time derivatives of the response variables, 
which are extracted to form the Inertia matrix, all load distributions such as 
Eqs. (4l) and (44) are calculated and used directly In the spanvd.se Integra- 
tions required for the various response excitations. Thus, the major necessary 

complexity Inherent in force-integration, that of having to calculate expli- 
citly the load distributions. Is for the most part already satisfied. The 
remaining requirement for inrplementation is tp approximate the doubly differ- 
entiated responses (vibratory accelerations, ^ ) and include them in the load 
distributions . 

The general time-history solution, as Indicated by Eq. (26), implies that 
at axiy one time step when the excitations are being calculated, the accelera- 
tion is not yet known: The (RHS) of Eq. (26), which nominally excludes all 

^ *s, is first evaluated to calcvilate the ^ 's given on the (LHS) of this 
equation. Thus, to Implement the force-integration calculation, approxima-r 
tlons to the 's implicitly appearing In the (RHS) must be made. Fortunately, 
the force-integration method is generally "forgiving" of such approximations 
since loads due to vibratory acceleration usually act as dynamic corrections 
to the "pseudo-static" loads. A discussion and formulation of the extrapola- 
tion formula required to approximate these accelerations is contained in 
Appendix Ij briefly, for any response with an inherent natural frequency, <u, 
the acceleration can be approximated as: 


«« 


«« 




(131) 


where Is the Integration step size. With this approximation formula, Eqs. 

(4la, b, and c) 8uid (44a, b, and c) can be evaluated and the following expres- 
sions for the Internal concentrated blade loads outboard and adjacent to the 
juncture can then be calculated. In these expressions, unless otherwise indi- 
cated, the load distributions are understood to contain both Inertia and 
aerodynamic contributions: 
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(132a) 

B fj 

(132b) 


( 132 c) 

® 'll 

(I 32 d) 


(I32e) 

“'5, ' [’’Ozs'*' P)'5*'‘'j’ “ '’0x5*''5”y5j’] 

(I32f) 


Additionally, the push-rod force, Ppj^, appearing in Eq. (112) must 
similarly be approoclmated. Upon using the principal (linear) component of the 
flexbeam torsional stiffness, the push-rod force can be approximated by: 

PR 

Pz5'y8-y5p-»iR<'-^ll» ( 133 ) 

-PygU5-Z5^-z'5^(r-rR))]dr} 


where the subscript (r) denotes conditions at the root of the torslonally 
active portion of the blade. The above approximations together with the re- 
sults of the two preceding subsections are sufficient to explicitly formulate 
the details of the rediuidant analysis solution which is accoapllshed in the 
following subsection. 
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Mathematical Formulation of the Structxiral Redundancy 


The results of the previous subsections can be brought together to form 
a solution for all internal loads in the vicinity of the juncture. From the 
deflections of the Juncture and inboard torque tube attachment point (snubber) 
all concentrated Juncture loads in the torque tube except torsion moment are 
evaloiated by means of Eqs. (112) and (113) • All loads outboard of the juncture 
are known from interaction of the spanwise loadings and are calculable vising 
Eqs. (I32a through f). Consequently, by equilibrating the internal Juncture 
loads and imposing the consistency constraint on the elastic torsion deflec- 
tion of flexbeam , the flexbeam and thence, torque tube torsion moments are 
determined. First, however, some attention must be paid to the internet, bend- 
ing moments in the torque tube at the Juncture. These moments as calculated 
using Eq. (112) are oriented in the torque tube section plane at the Juncture 
and must be rotated back to the "5" coordinate system using the following co- 
ordinate transformation pair: 










(I3l*a) 


(134b) 


The primary purpose of this redmdant analysis is to obtain the incre- 
mental torsion moment exerted by the flexbeam to the torque tube, which 

is given by: 


AM 




-(M, 


»FB 




(135) 


But the flexbeam loads are related to the outer blade and torque tube 
loads by the equilibration of loads: 



(136) 




Use -of the coordinate transformation^ Eq,. (134a) , on the outer blade 
moments Eqs. (I32d through f), the torque tube loads, Eqs. (112 and 113) and 
the equlbratlon of loads given by Eq. (136) sufficiently determines all loads 
In the ”7" coordinate system except which can be determined, however, 

fj*om the torsion stiffness equation of the flexbeam . Upon Inclusion of the 
effects of flexbeam bullt-ln twist, equation can be rewritten as: 

AOj * 6j - ^Z|fb 

* “y7FB> “^7fb’ 


and solving for yields the final required result: 


^SxSx5pg+ ®SySv5pg+ ®*2^^5 fb 


( 138 ) 


This Incremental moment Is then Included as a concentrated torsion moment 
to the torsion response, Eq. (30) as is therein provided, to complete the 
redundant analysis. 


The section to follow develops the eqviations needed for aeroelastic eigen- 
solutions using a fixed azimuth approach. While the linearized equations derived 
do contain the bending redundancy, they omit the (essentially nonllnecur) torsion 
redundancy described In this subsection. 
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LINEARIZED FORM OP EQUATIONS 


For time-history solutions of the complete nonlinear equations the conrplete. 
explicit expansions of the inertial and aerod^amlc load distributions are not 
required. Ih fact, the li^llclt descriptions of these load distributions (i^s. 
(4l), (44), and (59)) actiially facilitate the formulation of these equations 
and eliminate the requirement for assessing orders of magnitude in order to 
achieve tractibility through simplification. On the other hand, the unique suad 
desirable featTires of elgensolutlons (l.e., relatively short computation time 
for solution, identification of all coupled mode frequency and damping charac- 
teristics and the availability of established analysis techniques for assessing 
the behavior of linear systems) are sufficiently attractive to warrant equation 
linearization using explicit expansions of these loadings. The purpose of this 
subsection, therefore, is to present the main results of this linearization. 

Ihe relative orders of magnitude of the various dynamic and elastic quantities 
are presented in Table I; ftom this tabxilation the various coefficients in the 

TABLE I 

Assumed Relative Orders of Magnitude of Quantities Appearing in Aeroelastic 
Dynamic Equations. 
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resulting escpanded equations of motion vere evalmted and retained only if they 
did not exceed As a result, some terms so retained remain nonlinear, but 

are generally of a quadratic form (i.e., *3a*%» etc.). However, 

these terms can be locally linearized about predeflected positions. Ihe ex- 
panded equations presented in Table I are consistently expanded to in the 
coefficients and have the quadratic nonlinear dynamic and elastic terms segre- . 
gated for clarity. 

Heasonable explicit expansions and linearizations of the aerodynamic load 
distributions, contained only Inpllcitly in the above equations, can be achieved 
within the context of various optional simplifying assumptions. The aerodynamic 
linearization actually developed for the analysis parallels the aerodynamic de- 
scription given in an above subsection so that the total linearization achieved 
should represent a true "local linearization" of the toteO. nonlinear equations 
at any point in time. 

Using Eq. (1.59) as a starting point the following terms eure seen to con- 
tain all the independent response variables whose perturbations must be made: 


Up = 7p- VySinQ-v^cos® 


(139a) 


Ut = 7V+ VyCOS®-%sin® (I39b) 

«qs + + ton ' (-^) (139c) 

M = (139d) 

where : 

x(l--l^) -/SU*-r4 (lUOa) 

= e + r^l - + Sj + /A 1^1 - -^jsin^ + 8cost|>j (lUOb) 
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+ ©[(y,-AV,,)q,^+AViq^i + y,o^] + Ux[(y;.- Aw|^^)q^. 




(lllOd) 



4 


and where U^j. is given by Eq. (50). 

Assuming that the partial derivatives of the airfoil aerodynamic coeffi- 
cients, c^, C(j, and Cjuc/U’ wi'th respect to angle-of-attack and Mach number are 
available, the following ejjpressions can be written for perturbational aero- 
dynamic load distributions : 
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(142) 






'Ax. 2 mo 


+ ^ MUpj 8Up 

'• * 2Cp.Ut--^Up+^“Ut1sUt 

- '<A,5USe + y,o^^^8[coseppj^- $inep»,jr 


(143) 


where : 


8Up = Svp- (VyCOS© -V 2 Sin®)S© -sin@8v^-cos®8^ 


(l44a) 
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8Ut = Sv-p- ( Vy sin0+‘^cos0)S® +cos©SVy- sinQS^z 


and where SVp, 6'^_> and are obtained from appropriate differentiations 

of Eqs, (l40; and (62b; or ( 69 ). 

Further explicit expansions of the perturbational airload distributions 
beyond this point would rapidly become overly tedious, and voluminous and would, 
thus, not serve much additional purpose. The above development together with 
appropriate straightforward perturbations of the component air velocities are 
sufficient to form the desired perturbational airloads. The perturbational 
inertia loads are obtained in a similar manner using Eqs. (4la, b and c) and 
Eq.s. (44a, b and c). 
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APEE3SDIX I 


QUADRATURE FORMULAE USED IN PROGRAM G^lOO 


The time-history solution of the fully nonlinear set of dynamic eqmtions 
of motion, Eq. (26), requires an appropriate set of numerical integration, 
extrapolation and differentiation formulae. Numerical integration in (nondi- 
mensional) time is required for the basic solution of the nonlinear equation 
set in a conventional step-by-step fashion. Furthermore, when various of the 
elastic modal response vsuriables are characterized by large natural frequencies 
12P), numerically stable integrations would require prohibitively small 
integration step sizes. For such cases it becomes reasonable to assume that 
the inertia (^) term is negligible compared with the elastic (o^q) term. 

This assumption permits eliminating the double integration of the inertia term 
and substituting a "quasi- static" solution which involves an appropriate 
extrapolation of the response variable. Finally, whereas Eq. (26) implies ^ ^ 
that the right-hand side is a function devoid of explicit linear terms in q 
(the highest derivatives of q), various explicit or implicit nonlinear terms 
in these derivatives may still exist. Two particular sources of nonlineeor 
double time derivative dependency are the estimation of dynamic loads for use 
in the redundant analysis, and the calculation of the second time derivative 
of angle of attack as required for the unsteady airloads formulation (see Refs. 
9 and 10 for a detailed description of this formulation). Because of such 
nonlinearity, the usual time-hlstoiry algorithm of solving for at a given 
time Tising knowledge of the lower derivatives at that same instant is disrupted. 
Hence, approximation to based upon a reasonable extrapolation of past 
values, is required for evaluation of these nonlinear tearms on the right hand 
side of Eq. (26). The following subsections describe the numerical algorithms 
used to satisfy these three solution requirements. 

Numercial Integration 

Double integration of the response variable accelerations, to obtain 
q and q at the next time step is accomplished using a variant of the Adams 
method (without correctors), m this variant the accelerations are locally 
assumed to be simple harmonics of their respective characteristic frequencies 
in order to integrate over each subsequent time step. Denoting this frequency 
as an "integration" frequency, (H, the following algorithms result: 


k+i 


« 

= Ql 


C0SaiA)|^ - COS25JA^ 
(D Slnte>Ai|/ 




/ ^ *” COScjAi/^ 

\ usinoJAi^ / ^k-i 


(I.l) 
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*»^ere Is the Integration step size. For each elastic (modal) degree of 
flreedom the integration frequency is taken to he its Inputted uncoupled natvtral 
frequency. Ilie usual pendular frequencies of a rotating, hinged, rigid beam 
are taken as the integration frequencies of the flapping and lead-lag degrees 
of freedom. 


Quasi-Static Solution 


03ie static solution on any one of the response variables, q_, is achieved 
by manipulating the simultaneous equation solution of Sq. (31») 


*• XI * • \ 

“Isi, • '."’ll' Ok- ♦k* 


(1.3) 


into the following approximate foim; 


where, in general and are not equal. Hie newly formed quantity, 
is used only to effect an extrapolation to the subsequent time step: 


# I , A A . 

Ok., ■ IZ* <0k- 0k-2> 


(1.5) 


q 


k+i 




( 1 . 6 ) 


The results of Eqs. (l.^) and (l.6) are then used in exactly the same way as 
those of Eqs. (l.l) and (l.2) to evaluate the right-hand (excitation) side of 
the equations of motion, Eq. (26). Furthermore, it shotild be stressed that 
these equations give valid approximations only for, and should be limited to, 
those cases where the integration frequency is large. 
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Extrapolation of q_ 


As -with the numerical integration scheme presented in an above subsection, 
extrapolation of ^ is achieved by assuming the accelerations to.be locally 
simple harmonic. The resulting extrapolation formula becomes: 


IMI 

» 2(cpscA»/^)qi^_,- qi5_ 


( 1 . 7 ) 
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APEEINDIX II 


PROGRAM g400 INPUT DESCRIPTION 

The required input to the program consists of the following maj'or punched 
card data blocks in order of loading: 

I. Airfoil Data 

II. Inertia, Elastic, Geometric and other Operational Data 
(in Loader Format) 

III. Blade Mode Shape Data 

IV. Harmonics of Variable Inflow 

Details for preparing the data for each of these blocks are given in the 
sections which follow. An additional section of this appendix provides infor- 
mation for facilitating program operation and Improving its efficiency. 


I. Airfoil Data 

This data block consists of tables of two-dimensional lift, drag and 
pitching moment coefficients each versus angle of attack for vario\as Mach 
numbers. Additionally, if unsteady aerodynamics are used, the static stall 
angles and linear coefficient slopes for both lift and pitching moment are in- 
cluded in this table. Provision is made for Inputting eind using only one set 
of aerodynamic coefficients; hence, it is assumed that the same airfoil section 
is used over the entire blade span. I^rovlslon is also made in the program for 
optionally using an analytic representation of the NACA 0012 airfoil; if this 
option is invoked it is not necessary to provide any airfoil data to the pro- 
gram. For usage of the analytic NACA 0012 airfoil option, however, the 
required input for this block of data must be a single card with blank or 
zeroed columns 1 and 2. For the general case, however, the airfoil data are 
loaded in three subblocks (corresponding to c^j, c^, and Cmc/4 data) using the 
following general format; note that the required punch format is indicated by 
the FORTRAN format infonnation in parenthesis: 


card #1: 




TITLE (optional) (I2,AT8) 
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card ^ 


cards 




J N M A(l) CL(1) A(2) CL(2) ...A(4) CL(4) 

(12,1057.0) 


^(5) CL(5) ... 

(F9.o,9F7.o) 


^..A(N) CL(N) ALSTAL DCLDAO 

(F9.0,9F?.0) 


wherej NZ is the ntunber of Mach ntunbers for which groups of c^ data are to be 
read inj TITLE is any (optional) identifying information. J is the number of 
data entries to be inputted into each such group. N is the nxamber of angle- 
of-attack -c. abscissa-ordinate pairs to be inputted; N is restricted to a 
maximum of 3^^ without, and 33 with the use of unsteady aerodynamics. M is 
the Mach number appropriate to the data group. A(i) are the N smgle- of -attack 
abscissae in degrees and CL(i) are the N lift coefficient ordinates. AISTAL 
and DCLMQ are, respectively, the static stall angle, in degrees, and the 
lift curve slope at zero angle-of-attack, in per degree; these items ar« needed 
only if the \msteady airloads option is invoked. 

Cards 2 and 2+ are repeated for each successively higher Mach number. A 
meocimum of 12 Ifech nrunbers is allowed and the lowest and highest Mach numbers 
need not define the total working range as the search technique uses the 
boxmdary data for tfe,ch numbers beyond the inputted range. Th\is, repeated data 
for zero and supersonic Mach numbers are not needed. The lowest Mach number 
inputted must contain an angle-of-attack range of from -l 80 ° to l 80 ° or from 
0 to l80° depending on whether or not unsymmetric airfoil data is being 
inputted; all higher Ifech number data need extend only from -30° to 30° or 
from 0° to 30 ° in a similar manner. 
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!Hie general format described above is repeated for the and cmQ/h 
subblocks in that order. The static stall angles and aerodynamic coefficient 
cvirve slopes at zero angle-of-attack are deleted for the c^ subblock. 


II. Inertia, Elastic, Geometric and Other Operational Data 

This data block includes those items used to define the more detailed 
dynamic features and/or those which are npst likely to vary from case, to 
cate. The format for these data is as follows: 


NN L DAIA.(L) nATA(Lfl)...DATA.(Lf4) (l2,lU,5P12.0) 


where: UN is the card word count, i.e., the number of data items on the card 

to be inputted, columns 1 and 2; NK must not exceed 5. L is the location or 
identifying number of the first data item on the card coltanns 3-6 right 
adjusted. DA,TA.(Lfi) represents the various data items on the card, columns 
7-18, 19-30, 31-42, 43-54, and 55-66, in floating point format. The locations 
or identifying numbers for the various data auid control items eire listed below 
along with definitions and other pertinent comments; note that data locations 
not assigned data are implied to be intentionally left blank and/or not inputted: 


Location Item 

1 flR 

2 R 

3 P 

4 

5 b 

6 e 

7 B 


Description 

Rotor tip speed, ft/sec. 

Rotor radius, ft. 

Air density, Ib-sec^/ft^. 

Speed of sound, ft/sec. 

Number of blades. 

Nondimenslonal offset distance of start 
of deformable and/or deflectable portion 
of rotor blade, e/R, 

Tip loss, used to define equiveklent mo- 
mentm area and three-dimensionality 
corrections to computed two-dimensional 
airloads near the blade tip. 
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Location 

8 

9 

10 


11 

12 

13 

14 

15 

16 
17 


Item 

Description 

NSEG 

Number of blade segments used to define 
spamd.se variable arrays. 

Lif 

Azimuth Increment used in the numericeJ. 
integration of the dynamic equations, 
deg. (See section on general informa- 
tion for efficient program usage. ) 

% 

Number of "flap trials", l.e., maximum 
number of rotor revolutions for which 
the blade time-history will be computed 
in an attempt to obtain convergence to 
periodicity. If a transient response is 
desired for only a portion of one rotor 
revolution the program will compute a 
time-history solution for any nonzeiro 
fractional Np value inputted. An Iden- 
tlceOJLy zero value will cause the time- 
history solution to be by-i>assed entirely. 


Flapping tolerance to within which the 
aeroelastic/dynamic responses must repeat 
on successive revolutions in order for 
the motion to be considered converged to 
periodicity. The tolerance applied to 
lead-lag motion is eqttal to ?Cp. 

At 

Longitudinal cyclic pitch, coefficient 
of minus cos^ term in Fourier expansion 
of blade control pitch angle, deg. 

% 

-■-s 

Lateral cyclic pitch, coefficient of 
mlnois sln^ term in Fourier e:^ansion of 
blade control pitch angle, deg. 

®.75R 

Blade collective pitch angle as defined 
at the 75^ radius, deg. 

\ 

Mean rotor Inflow ratio. 

V 

Forward flight velocity, kts. 

a 

Rotor solidity (■ bc/nR). 
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Location 


Item 


18 


19 


20 


21 

22 


23 


24 


25 




(Control) 


N 


cut-out 



cut-out 


01 


^ 'Sprint 


e 


6 


Description 


Increment added to all values of 
obtained from tabulated airfoil data Or 
from the analytic MCk 0012 data. Air- 
foil data generally correspond -to smooth 
wind tunnel models and is often 
used to adjiist for the higher drag of 
production blades; a commonly used value 
of AccLq is 0.002. 

Make greater than zero (1. ) for first 
case or when new blade modal data are 
to be inputted. Program automatically 
sets this control number to (-1.) after 
each loading of modal data. 

Number of blade segments, starting at 
inboard end and defining the cut-out 
region, for which the lift and moment 
coefficients are set to zero. 


Ihe drag coefficient used on the first 


N 


cut-out 


segments . 


Built-in linear blade twist angle; i.e., 
difference between tip and root built-in 
angles, positive when tip angle is 
greater (L.E. up) than root angle, deg. 


Azimuth increment used to present printed 
output of various pertinent aerodynamic, 
dynamic and elastic load distributions 
as -well as aeroelastic responses and 
stresses, deg. Biis input quantity 
should be an integral multiple of loca- 
tion 9 j deg. 


"Direct” value of pitch-flap coupling 
(a A0/dP)» -A- nonzero value will sup- 

press a calculation of this quantity 
from the inputted pitch-hom/push-rod 
geometry. 

Pitch-lag coi;^ling (» Ao/as)* 
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Location 

26 

27 

28 . 

29 

30 

31-34 

37-39 

Uo-h2 

43 

44 


Item 

Description 

^LD 

Viscous lag 

sec. 


damper coefficient, ft-lb- 


c 

NFM 

NEM 

NTM 


“Vi 







Blade chord if chord is constant, 
otherwise omit, ft. 

Number of flatwise bending modes to be 
used (4 max). 

Number of edgewise bending modes to be 
used (3 max). 

Number of elastic torsion modes to be 
used (3 meoc) (note that the total num- 
ber of dynamic degrees-of- freedom is 
limited to 10. Ihiis, if the articula- 
tion degrees-of- freedom, p and/or 6, 
axe used, or if the rigid feathering 
d.o.f. is invoked, NFM, NEM and N!M will 
be automatically changed to keep the 
total degrees-of-freedom to no more than 
10.). 

Flatwise modal frequencies, nondlmen- 
sional with respect to fi, in- ascending 
modal order. 

Edgewise modal frequencies, nondlmen- 
slonal with respect to 0, in ascending 
modal order. 

Torsion modal frequencies, nondlmen- 
slonal with respect to 0, In ascending 
modal order. 

Second harmonic cyclic pitch coefficient 
of minus cos in Fourier e:^)ansion of 
blade control pitch, deg. 

Second harmonic cyclic pitch coefficient 
of minus sin 2«/f in Fourier expeuision of 
blade control pitch, deg. 
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Item 

g 

(Control) 

(Control) 

(Control) 

^root 


(Control) 


Description 

p 

Acceleration due to gravity, ft/sec j a 
negative value implies inverted flight. 

Make nonzero (1. ) if airfoil data for a 
nonsvmmetric airfoil are to be used. 

Make nonzero (1. ) if the total ( transi- 
ent ) time-history is to be outputted; 
i.e., responses calculated before con- 
vergence to periodicity is obtained. 

Make nonzero (1, ) if the modal responses 
and hub shears and moments axe to be 
(negative) Fourier anetlyzed after peri- 
odicity has been obtained. 

Torsional rigidity at the blade root, 
Ib/ft^. 

Torsional spring connecting root of blade 
to fixed structiire to represent control 
system flexibility, ft-lb/rad. A non- 
zero value will automatically introduce 
the rigid-body feathering degree-of- 
freedom as an addition "torsion mode". 
Note that this "mode" will inherently 
couple with the NTM normal elastic tor- 
sion modes at frequencies both below and 
above the inputted torsion frequencies 
(locations kO-k2); hence, a smaller 
integration inte 2 nral will generally 
be required. 

Make nonzero (1. ) to simplifir the nvmer- 
ical spanwlse integration tecliniques 
from the nominal trapezoidal rule to 
rectangular (Eulerlan) form. Rectangu- 
lar integration is effected by setting 
the quadrature nmbers \ised for spanwlse 
integration equal to the Inputted seg- 
ment lengths (loc. (100-114)). Usage of 
this option is recommended for blade 



Location 


52 


53 


54 

55 


56 

57 

58 


59 


Item 


(Control) 


(Control) 


(Control) 


a 


^req 


PP. 


req 


®Lift 


*PP 


Description 

configurations with slgnlflcEint 
discontintdties in spanwise properties 
(e.g., counterweights, tip weight, step 
twists, etc.)* 

Make nonzero (1. ) to output modal 
Integration constants used in the 
elgensolutlon and, to a limited extent, 
in the time-history solution. 

Ifeke nonzero (1. ) to load prepunched 
(vorticity) induced velocity distribu- 
tions . 

Make nonzero (1. ) to use the Induced 
velocities loaded as per location 53* 

Shaft angle of attack, deg. Ihis input 
item serves a dual role; when variable 
inflow is used (either vorticity induced 
or Glauert momentum) is used to define 
Xpgm* Also, when a major Iteration is 
to be performed (nonzero location 60 ) 
and when trims on propulsive force are 
deactivated (zero location 59) "the pro- 
gram will trim to this inputted shaft 
angle. 

Requested value of lift to be used in 
major iteration, lb. 

Requested value of propulsive force to 
be used in major iteration, lb. 

Tolerance on lift for major iteration, 
lb. A zero value deactivates trimming 
on lift. 

Tolereuice on propulsive force for major 
iteration, lb. A zero value deactivates 
trimming on propulsive force. !l3ie auto- 
matic trim calculation (major Iteration) 
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Description 


must trim either to a required propulsive 
force or to a required shaft angle of 
attack; therefore a deactivation of 
trim to propulsive force automatically 
directs the trim calculation to trim to 
requested shaft angle, location 55* 

Maximum number of major iterations to 
be made in an attempt to achieve trim. 

A zero value will deactivate the major 
iteration; a negative value will acti- 
vate a stall avoidance calculation if, 
when attempting to trim, the rotor be- 
comes stalled. (See section on general 
information for efficient program usage. ) 

Built-in precone angle, deg. An identi- 
cally zero value denotes a rotor blade 
articulated in flapping; similarly, a 
finite nonzero value signifies a blade 
nonarticulated in flapping. 

Requested value of pitching moment for 
major iteration, Ib-ft (positive nose 
up). 

Tolerance on pitching moment for major 
iteration, Ib-ft. A zero value deacti- 
vates trimming on pitching moment. 

Requested value of rolling moment for 
major iteration, Ib-ft (positive port 
side up). 

Tolerance on rolling moment for major 
iteration Ib-ft. A zero value deactives 
trimming on rolling moment. 

Initial condition on articulated flap- 
ping angle, rad. 

Initial condition on (nondimensional) 
articulated flap angle rate. 



Item 

6 

* 

6 

(Control) 

(Control) 

(Control) 

(Control) 

(FOR) 


Description 


Initial condition on articulated lead- 
leig angle, rad. (positive forward). 

Initial condition on (nondimensional) 
articulated lead-leig rate. 

Built-in prelead angle, deg. An iden- 
tically zero value denotes a rotor 
blade articulated in lead- lad; similarly 
a finite nonzero value signifies a blade 
nonarticiilated in lead- lag. 

Make nonzero (1. ) to harmonically 
analyze and output harmonics of flat- 
wise stresses. 


Make nonzero (1. ) to harmonically 
analyze and output harmonics of edgewise 
stresses. 


Jfeke nonzero (1. ) to harmonically 
analyze and output harmonics of torsional 
stresses. 


Make nonzero (1. ) to bypass the eigen- 
solution. 


Initial condition on rotor azimuth, deg. 
This value, like all inputted initial 
conditions, is likewise used to define 
the loadings on the blade about which 
linear perturbations are to be taken in 
the eigensolution. 

Factor in momentum inflow equations to 
account for dual (coaxial) rotor opera- 
tion: 




Note that C^ in the above equation is 



Location 


Description 


79 


80-82 


83 


Item 


(Control) 



SR 


the usual thrust coefficient as defined 
for a single rotor. An identically zero 
value of (PCR) defaults to a value of 
2.; a value between 1. and 2. will de- 
note a coaxial rotor operating some- 
where between hover and high speed 
forward flight. 

Generalized Glauert (momentum derived) 
variable inflow option. A zero value 
deactivates usage. A value of 1. 
causes the inputted induced velocity 
components to be used, as inputted; a 
value of 2. causes the inputted values 
to be used initially, and then varied 
in trim calculations in place of control 
angles; a value of 3. causes the control 
and shaft angles to be fixed and the 
induced velocity components to be varied 
only to satisfy momentum equations in 
any requested major iteration. If the 
value is 1. , the induced velocity com- 
ponents will be varied to satisfy mo- 
mentum considerations in addition to 
the usual trim calcxilation. 

Initial conditions on the "momentum" 
induced velocity components comprising 
a Glauert-like variable inflow descrip- 
tion. Note that the "vorticity" vari- 
able inflow (controlled by locations 53 
and 5^) and the momentum variable inflow 
can be used separately or simvLLtaneously. 

Sample rate for Transient Spectral 
Stability Analysis (TSSA). (See Appen- 
dix III for a discussion of this tech- 
nique.) Every (SR)'th point in a 
transient time-history is saved for use 
in a TSSA. A zero value bypasses the 
TSSA. 
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1 



Location 

Item 

Description 

84-86 

(Control) 

Channel selection for each of three 
available for the TSSA. Bie channels 
available are; 1-4, flatwise bending 
modal responses; 5-7> edgewise bending 
modal responses; 8-10 torsion modal 
responses; 11, articulated rotor flap- 
ping; 12, articulated rotor lagging; 

13 > l4, and 15, blade tip vertical, in- 
plane and torsion deflections, respec- 
tively. 

87 

“l 

Lower bo\md of frequency band chosen for 
TSSA nondimens lonal with respect- to Q . 

88 


Upper bound of frequency band chosen for 
TSSA nondlmenslonal with respect to 0. 

89 

(Control) 

Initial estimate of percentage of total 
transient data used in each time dls- 
placed data sajsple block in TSSA. 

90 

(Control) 

Number of transient (time displaced) 
Fotirier coefficient calculations made 
to establish modal damping in TSSA; 
maximum value is 200. 

91 

FREQ 

Nianber of desired resonant frequencies 
to be extracted from frequency band 
defined by locations 87 and 88. 

92 

(Control) 

Resonant frequency identification 
criterion. Make (0. , 1. ) to identify 
resonant frequencies using criterion 
of (maxiiiium Fo\irler Transform amplitudes, 
maxlirnim F.T. anplltudes + deletion of 
"Image** frequencies). Note that image 
frequencies result from the Floquet-llke 
response characteristics of rotors in 
forward flight. 

94 

^vse 

Critical (viscous) damping ratio used to 
approximate structured, dashing in edge- 
wise bending modes. 
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Location 


Item 


Description 


97 


98 


99 


100-114 


115-129 

130-144 


®^equ 


(Control) 


(Control) 


AX 


% 

®Ba 


Linear equivalent blade (nonlinear) 
twist angle defined similar to location 
22 , deg. A nonzero value for 0 ]^ is 
required to use the inputted nonSSnear 
twist arrays and is used to calculate 
that portion of the aerodynamic inflow 
velocity at the 3/4 chord position. Up, 
accruing from radial flow and twist. 

Make nonzero (1. ) for stress calculations 
using the mode deflection method. Zero 
value defaults to force integration 
method. 

Location used to end a case or series 
of cases. Make (+1.) to end the Loader 
Format data block for the case defined 
by the Loader data and load additional 
cases at the conclusion of that case. 
Make (-1. ) to end the Loader data and 
read no further cases. In both insteinces 
the word count, LL, (see beginning of 
this section above) should be (-1). 

Note; this entry must appear singly on 
an input card, and that card must be the 
last card for the case. 


Nondimensional blade segment lengths, 
in order from root to tip, mav-iTmim of 
15 values, starting from the offset lo- 
cation. Accxiracy is generally improved 
if the last segment is small (s0.03). 

Mass of each blade segment, lb-sec^/ ft. 

Aerodynamic built-in nonlinear twist 
angle distribution, deg. Since collec- 
tive angle is defined at the 75^ span 
location, should have a zero value 
at 75'^ span. Should the structural 
twist angle distribution differ from 
9g^, their appropriate data must be 
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Location 


Item 


Description 


145-159 

160-174 

175-189 

190-204 

205-219 

220-234 

235-249 

250-264 


C 




K. 


t/m 



^lOc/4 


loaded into locations 690-704; othend.se, 
0D will be used for both aerodynamic 
am structural applications. Nonlinear 
twist distributions will be used only 
if a nonzero value is inputted into 
location 97, 

-*-equ 

Blade chord at center of each blade 
segment (use for nonconstant chord 
blades only), root to tip, ft. 

Constants relating (nondimensional) 
flatwise second derivative to flatwise 
stress (s» Ec/r)j,, evaluated at center 
of each segment, root to tip, psi. 

Constants relating (nondimens ioneJ. ) 
second derivative to edgewise stress 
(=Ec/R)g, root to tip, psi. 

Constants relating torsional moment to 
torsional stress, root to tip, in"^. 

Chordwise mass radii of gyration of 
blade segments about elastic (reference) 
eixis, root to tip, nondimensional with 
respect to R. 

Thicknesswise mass radii of gyration of 
blade segments about axis perpendicular 
to chord line and throxjgh the reference 
axis, root to tip, nondimensional with 
respect to R. 

Area radii of gyration about elastic 
eixis, root to tip, nond imensioneil with 
respect to R. 

Distances from elastic axis forward to 
airfoil quarter chord position, root to 
tip, nondimensional with respect to R. 


107 



Location 

Item 

Description 

265-279 


Distance from elastic axis forward to 
airfoil section mass centers, root to 
tip, nondimenslonal with respect to R. 

280-283 


Initial conditions on i’th flatwise 
bending mode deflections. 

286-289 

* 

GWj^ 

Initial conditions on i’th flatwise 
bending mode (nondimenslonal) rates. 

292-294 I 


Initial conditions on k’th edgewise 
bending mode deflections. 

295-297 

* 

Initial conditions on k'th edgewise 
bending mode (nondimenslonal) rates. 

298-300 


Initial conditions on j’th torsion mode 
deflections. 

301-303 

* 

Initial conditions on j’th torsion mode 
(nondimenslonal) rates. 

600-6IU 


Section modiilil for flatwise bending, 
root to tip, in^. Note that the pro- 
ducts of (l/c)p and Kp (locations I6O- 
174 ) must equal the bending stiffness 
distribution, EI/r. 

615-629 


Section modulll for edgewise bending, 
root to tip, in 3 . See remarks above 
for flatwise bending section moduli! , 

630-639 

(Control) 

Stress selection nimiber for each of (at 


mcudjnum) 10 ch€umels whose time-history 
data strings are saved and used for 
automatic plotting. Value is determined 
by formula; J x 100 + N, where J a 
( 0 . 1 , 2 , 3 ) as stress is (flatwise, edge- 
wise, torsion, flexbeam torsion) and 
N ■ segment nxmiber. (input locations 
provided, but option inoperative.) 
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Location 


Item 


Description 


640-649 

660-674 

675-689 

690-704 

705-719 

720-734 

735-744 


STRSCALE 




EB^ 




40b 



Vertical axis scaling for the automatic 
plotting of each stress channel selected 
in locations 630-639* psi x ^0^/±a.. 
(input locations provided, but option 
inoperative.) 

TorsioneJ. stiffness (to be multiplied 
by twist rate squared), as defined in 
reference 4, Ib-ft^. 

Torsion to edgewise elastic coupling 
stiffness (to be multiplied by twist 
rate), as defined in reference 4, Ib-ft^. 

Structureil built-in nonlinear twist 
angle distribution, root to tip if 
different from aerodynamic twist, deg. 
See remarks above for aerodynamic 
built-in twist, location 130-144. 

Built-in (struct'ural) twist angle 
change per segment length distribution, 
root to tip, deg. Note that this item 
is a direct statement of the b\iilt-in 
twist rate distribution, e^; if al l 
values of this distribution are inputted 
as zero, the twist rate distribution is 
computed internally using nmerical 
methods from the inputted twist angle 
distributions, locations 130-144 or 
690 - 704 , as appropriate. 

Distances from reference (elastic) axis 
forward to edgewise bending neutral 
axis, root to tip, nondimensional 
with respect to R. 

Elexbeam plate bending stiffness dis- 
tribution, Ib-ft. Note that all items 
in locations 735 through 824 are re- 
quired only for a redundant anedysis of 
the torque-tube ( flexbeam span of bear- 
ingless rotor (CBR.) configurations) (see 
location 991 )• 
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Location 


Item Description 


7 ^* 5-754 

GJpB 

Flexbeam torsion (St. Venant) stiffness 
distribution, Ib-ft^. 

755-764 


Flexbeam section width (chord) distri- 
bution, ft. 

765-774 


Torque tube flatwise bending stiffness 
distribution, Ib-ft^. 

775-784 


Torque tube edgewise bending stiffness 
distribution, lb. ft^. 

785-794 

%Ti 

Mass of each torque tube segment whose 
span is defined by location 100-109, Ib- 
sec^/ft. 

795-804 


Distances from reference axis forward to 
torque tube mass centers whose masses are 
given in locations 785 - 794 , nondimens ional 
with respect to R. (input locations 
provided, but quantity is not used.) 

805-814 


Section modulii for torque tube flat- 
wise bending, in^. (input locations 
provided, but quantity not used. ) 

815-824 


Section modxilii for torque tube edgewise 
bending, in^. (input locations provided 
but quantity not used. ) 

957 

h 

Height of positively thrusting rotor 
from groiond or wind t\mnel floor for 
purposes of evaluating Heyson correc- 
tions to rotor angle of attack (see 
Ref. 11 ), ft. 

958 

%T 

Height of wind tunnel test section for 


purposes of evaluating Heyson correc- 
tions, ft. Note that zero values for 
the wind tunnel test section dimensions 
Implies that ground effect corrections, 
rather than wind tunnel wall corrections, 
are to be made. 


no 



Item 


%T 


CASE 



KDEL3I 


NDEL30 


■pR 


yiOpj^ 


Description 


Width of wind ttuinel test section for 
purposes of evaluating Heyson correc- 
tions, ft. 

Case number. 


Built-in elevation angle of push-rod, 
attachment point, measurable when 0 
0, deg. 


s75 


Segment number of inboard attachment 
point of pitch input structural members 
(Pitch horn, feathering ciiff or torque 
tube) to blade spar. 

Segment number of outboard attachment 
point of pitch input structural member 
to blade spar. Note that if this struc- 
tural element is attached to the blade 
spar at a single point (cantilevered 
configuration) NDEL3I and NDEL30 must 
both be inputted with the same appropri- 
ate value. For bearingless rotor (CBR) 
applications NDEL30 serves the addi- 
tional function of defining the outer 
limit of the flexbeam; this limit is 
taken as the inner boTuidary of the 
RDELSO’th segment. 

Radial location (blade station) of pus.x- 
rod to pitch input structural member 
attachment point, in. 


Location forward of feathering axis of 
push- rod to pitch input structural mem- 
ber attachment point, in. Rote that 
input items 977-980 together with the 
inputted modal data provide the geomet- 
ric data from which pitch- flap, pitch- 
flat and pitch- edge coupling are calcu- 
lated. These calculations are bypassed 
if either of locations 978 or 98 O are 
inputted identically zero. 



Item Description 

Nojig Segment nmber of inneimost end of the 

feathering flexure; defaults to 1. 

9n g Built-in twist angle of the flexheam, 

positive leading edge up at outboard 
end of fl^beam, deg. Note that 83^^ 
is defined relative to the inputted 
twist angle distribution, which in turn 
is defined for zero collective angle. 


GJ. 


FB 


Average torsional stiffness of the flex- 
beam Ib-ft^. 




(Control) 


(Control) 



(Control) 


Equivalent or critical torsion stress 
coefficient for flexbeam ; i.e., stress 
per unit average torsion strain as de- 
fined by total flexbeam twist, in-psi. 

ifeke nonzero (1.) to include the " wob- 
ble" mode for cantilever torque tube 
configured bearingless rotors. 

When finite snubber stiffnesses axe 
inputted (locations 98? and 988) the 
stiffnesses are_ assumed to be aligned 
(vertically and horizontally, torque 
tube flatwise and edgewise) as location 
986 is (0., 1.) 

Snubber vertical (torque tube flatwise) 
spring rate, Ib/ft. 

Snubber horizontal (torque tube edgewise) 
spring rate, Ib/ft. 

When the torque tube flexbeam red\uidant 
analysis is Invoked (location 991) the 
torque tube is assvuned to be in ( tension, 
compression ) as location (989) is 
( 0 ., 1 .). 



Location 


990 

991 


992 


993 


994 


995 


996 


Iton 

^Oijij 

(Control^ 


(Control) 




«W0B 


^VOB 


(Control) 


Description 

Segment number of Innermost end of 
torque tube; defaults to 1. 

Make nonzero (1. ) to activate the re* » 
dundant analysis option. When the 
redundant analysis is activated, the 
flexbeam and torque tube stiffness and 
mass Inputs, locations (735) through 
(794), must be appropriately nonzero. 

Used to weigh the effectiveness of non- 
linear |[^EI bending excitation of torsion 
Over the flexbeam span; the effective- 
ness is taken to be (100?^, 0^, 100^ (x) 
pseudo-torsion mode shape as location 
( 992 ) is (0., 1., 2.)). 

Flexbeam aspect ratio parameter 

g ... ... 

(a — Vl.5 (1-1/)) to detesrmine plate- 
like deflection for torsion "pseudo- 
mode”. A zero value gives a rod-like 
deflection shape. 

Spring rate of rotatlonail spring con- 
necting outbosird end of torque tube to 
the blade spar about chordwise axis, 
ft-lb/rad. A zero value implies zero 
stiffness for pinned-pinned torque tube 
configurations and infinite stiffness for 
snubbed configiorations . 

Spring rate of rotational spring connect- 
ing the outboard end of the torque tube 
to the blade speir about an axis perpen- 
dicular to the chord line, ft-lb/rad 
(see above). 

Analytic (static) airfoil option . Make 
nonzero (1. ) to use the built-in axialytie 
approximation to the static NACA 0012 
airfoil data. 
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Location 


997 


998 


999 


1000 


1001 


Item 

(Control) 


Ail(p389 


K , 4 . 

plot 


(Control) 


'^S.75 


Description 


Dfasteeidy airfoil data option . Make 
nonzeTO to use the synthesized 

unsteady airfoil data technique described 
in reference 10, This can be used with 
either the explicitly inputted airfoil 
data or the built-in analytic NACA 0012 
static airfoil data. For nonzero values 
greater than unity, the "cutoff Mach 
number" (Mach number above which the un- 
steady data is disregarded in favor of 
quasi-static data) is taken to be the 
inputted value minus one (1.) rather 
than the built-in default value of 0,6, 
However, only nonstandard cutoff Mach 
numbers between 0,1 and 0,95 are 
accepted. 

Azimuth increment \ised to generate 
punched card output to be used as input 
data to UlRC program F389 for con 5 >utlng 
variable inflow, This option is by- 
passed with a (O.) input value, A 
description of the format of these out- 
put piuiched cards is contained in 
reference 12, 

Scale for automatic plots of varioTis 
time-history data strings, see locations 
84-86 and 630-639* (input location pro- 
vided, but option is inoperative, ) This 
option is currently inoperative but 
would be bypassed with a (0.) input 
value, when operational. 

Make nonzero (1, ) to activate usage of 
the tabvilated time-histories of incre- 
mental control angles. 

Number of abscissa-ordinate point pairs 
used to define time-history of (t); 

calculation of this time-history ie by- 
passed with a (0,) value. 


Il4 



Location 


Item 


Description 


1002-1050 75 Table of A 0,75 abscissa-ordinate pedrs; 

U(A0,75) ■ deg.; U(t) ■ sec, 

1051 Number of abscissa-ordinate point pairs 

s used to define time-history of (t); 

zero value bypasses calculation. 

1052-1100 AAi Table of abscissa-ordinate point 

pairs; U(Ai ; • deg; U(t) » sec, 
s 

1101 Number of abscissa-ordinate point pairs 

s used to define time-history of ABig(t); 

zero value bypasses calculation. 


1102-1150 Table of abscissa-ordinate point 

® pairs; U(& ^ « deg; U(t) » sec. 



III. Blade Mode Shape Data 

Included in this data are the radial distributions of the blade 
(uncoupled) flatwise, edgewise and torsion normal mode shapes and their 
derivatives. Ihese quantities must generally be provided from an external 
source such as Uhlted Technologies Corporation program E 159 or an equivalent, 
in- the following punched card format : 
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irtiere: NFM, NEM, and NIM are, respectively, the numbers of flatwise bending, 

edgewise bending and torsion normal modes whose mode shapes and derivatives 
are to be inputted. NSEG is the number of blade spanwise stations for which 
the inputted modal data are defined. F(l) are each of the. below listed modal 
functions (defined at the i'th spanwise stations) and I must be either 1, 6 
or 11 reflecting the requirement of five entries per card until NSEG entries 
are. made for each F function inputted. Ike modal functions must be loaded 
in the following order: 




V, 


^2 


(NSEG values, root to tip) 

^ If If If If If ^ 


(NEM flatwise modes) 





(NEM edgewise modes) 


(NM torsion modes) 


IV, Variable Bannonic Inflow 

If location 53 of the Loader block of operational data is nonzero, the 
following block of variable hsunnonic inflow is loaded in: 

cetrd #1: 


card #2; 


card. #3; 


card #3 + (NHAEMl/5): 



Subsequent cards repeat the pattern wherein I, the spanwlse station index, 
varies from 1 to NSEG. NHAEML is the number of harmonics of inflow to be 
loaded, and XLAMO(l), XLAMC(l,n) and XLAMS(l,n) are, respectively, the zeroth, 
n*th cosine and n'th sine components of harmonic Inflow at the I'th radisil 
station, wherein a positive Fourier series is assumed. 13ie Inflow is defined 
positive up and has the units of ft/sec. 


V. Multiple Case Runs 

Ihe above described data set-up defines the correct ordering of required 
data blocks for a general case. When multiple cases are run the second and 
subsequent cases utilize most of the data Inputted for the first case. USie 
following rules apply to the running of multiple cases: 

1. Airfoil data is loaded only for the first case; ckU subsequent cases within 
the run Tise the same tabular data, if analytic data is not used. 


H7il 




2. Only those items within the operational (Loader) data which are to he 
changed from case to case need be inputted. 


3. Item 99 of the operational data controls the running of subsequent cases; 
a (+1.) value causes a subsequent case to be loaded whereas a (-1.) value 
terminates the computer run after the current case. 

4. Unless otherwise specified (by a +1. value for operational data item 19) 
the Inputted modal array data block is used for all cases within the run eind, 
hence, no subsequent input of this data need be made. 

5. Similarly, unless otherwise specified (by a +1. value for operational 
data item 53) the Inputted harmonic variable inflow data block is used for 
all cases within the run and, hence, no subsequent input of this data need 
be made. 

6. Cperational data items 19 and 53 discussed above are both automatically 
set to zero at the conclusion of the data input for every case. 

7. Terminal conditions on the blade azimuth angle, item 77 > and on the 
degrees-of- freedom, items 67-70, and 280-303, for any case are carried over 
as initial conditions on these quantities for the subsequent case. Thus, 
for some applications, e.g., investigations of unstable responses, it would 
be appropriate to reinitialize these items on the subsequent cases. 


General Information to Facilitate Operation of Erogram 
And Improve Efficiency 

Aside from considerations of the actual aeroelastic parameters describing 
the blade configuration, which are covered in the above sections, additional 
attention should be paid to the mechanics of obtaining efficient numerical 
solutions of the dynamic equations. In this regard, there arise two basic 
areas of concern wherein this section should be of assistance. The first of 
these areas is the proper selection of parameters for efficient temporal 
numerical integration of the dynamic equations (flapping or minor iterations) 
and the second is the proper selection of parameters for effecting a satis- 
factory rotor trim (major iteration). The following subsections provide 
information for making proper parameter selection in each of these areas. 

Temporal Numerical Integration - As is discussed in Appendix I temporal 
integration of the higher differentiated response variables to obtain the 
lower ones is achieved in the g 400 program using a variant of the Adams inte- 
gration algorithm (see equations I.l and 1.2), The selected algorithm is 
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defined by means of the azimuthal integration step size, ^\|f, and the 
integration frequency, 

Ihe integration step size should be eui integral divisor of 36O; a proper 
choice depends on the maximum coupled frequency inherent in the various 
aeroelastic responses. A reasonable upper limit for ^\(r is 45 divided by the 
maximum such frequency in per rev. Values of greater than this upper limit 
will compromise the integration accuracy and, for sufficiently large values, 
will cause the conputed responses to develop "n\mierical" instabilities. As 
a corollary, a check on any response which is predicted to be unstable by the 
analysis, is to rerun the case with a reduced integration step size to test 
for the possibilitjr of the unstable response being merely a nximerical insta- 
bility. 

For each response degree-of- freedom a different integration frequency, 

CD, is used in the integration algorithm; this frequency is, for each of the 
elastic modes, the respective inputted natural frequencies (locations 31-34, 
37-39> and 40-42). Ihe usual pendular frequencies of a rotating, hinged, 
rigid beam are taken as the integration frequencies of the flapping and 
lead-lag degrees-of- freedom. In addition to defining modal stiffnesses and 
integration frequencies, the inputted frequencies serve yet another purpose. 

As noted above, the proper value of integration step size, ^\(r, varies 
inversely with the maximum modal frequency. Hius, run times (caused by 
reduced step size) will significantly increase as any one modal frequency 
increases. Since any degree-of- freedom exhibiting a large natura^j^ ^ frequency 
tends to respond quasi-statically, i.e., as if the acceleration (q) term 
were negligible, a reasonable approximatlon^to the response calculation is 
to avoid the numerical integration of the q term entirely and treat the 
response quasi-statically. This option can be invoked for any such high- 
frequency mode by inputting the appropriate frequency negatively ; a negative 
sign will not affect the proper usage of the frequency in the calculation of 
the dynamic equations. Note that this optional response calctilation can be 
Invoked singly pr in combination for any of the elastic modal responses, 
(negative values in any of locations 31-34 and 37-42). 

Hub Force and Moment frim - Operation of the trim'or major iteration feature 
of the program is controlled by input locations 55 thru 60, 62 thru 65 and, 
for some applications, 79* She main control for the major iteration is loca- 
tion 60, the number of major iterations, Nj^. A zero value cavises the major 
iteration feature to be completely deactivated. On the basis of past usage, 
a reasonable range for this input appears to be from 5 to 10, depending on 
the tightness of the convergence tolerances selected and the "goodness'- of 
the initial guesses on the control parameters. Convergence of the major 
iteration is adversely affected by any lack of convergence of the responses 
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to periodicity vithin each trim iteration and by incursion of the rotor into 
a significantly stalled flight regime. Note that a negative value of Nj^ 
will activate a stall avoidance calculation wherein the controls will be 
perturbed to reachieve unstalled flight rather than to achieve the prescribed 
trim condition. Should a ma^jor iteration fail to converge within any one run 
the last used control angles and initial conditions on response variables axe 
generally available in output punched card form and/or partially in the out- 
put printed records of each major iteration (see Appendix III) for use in 
subsequent major iterations. 

The trim iteration is operationally flexible as to what hub loads it 
will drive to requested values. Generally, the various requested hub loads, 
lift, propulsive force, pitching and rolling moment are activated in turn 
by specifying nonzero values for each of their respective tolerances. More 
specifically, the following table describes the optional combinations of hub 
loads and rotor shaft angle available with the G^tOO trim capability: 


TABLE II 

SUMMARY OF OPTIONAL BASIC TRIM COMBINATIONS 


Option 

Uft 

Prop. Force 

Htch. Momti 

Roll Momt. 

®.75R 

AV 

-*-s 

% 

■‘■s 

Of 

s 

1 

S 

U 

U 

U 

V 

F 

F 

s,(v) 

2 

s 

S 

U 

U 

V 

F 

F 

V 

3 

u 

U 

S 

S 

F 

V 

V 

S,(V) 

4 

s 

U 

S 

S 

V 

V 

V 

S,(V) 

5 

s 

S 

S 

S 

V 

V 

V 

V 

6 

u 

U 

S 

U 

F 

V 

F 

S,(v) 

7 

s 

U 

S 

U 

V 

V 

F 

S,(V) 

8 

s 

S 

S 

U 

V 

V 

F 

V 


where : 


F: Control parameter kept fixed 

V: Control parameter varied 

S; Trim parameter specified £ind trimmed to 
U: Trim parameter unspecified and ignored 
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Note that for each option available, the Individual trim pareuneter Is 
activated by Inputting a positive, nonzero value for the respective parameter 
tolerance. For each option the trim calculation trims either to the specified 
propulsive force or the Inputted shaft angle, ffg (location 55) as the propul- 
sive force tolerance, epp. Is flid-te or zero, respectively. A reasonable 
choice of tolerance values for the hub force and moment loads can be obtained 
from considerations of the helicopters Inertia properties and acceptable de- 
viations of the load factor from, unity, and of the pitching and rolling accer- 
atlons from zero. 

Ihe trim calculation can be c^so used In conjunction with the Glauert 
variable Inflow. Location 79 controls the usage of this slnqpllf led form of 
variable Inflow. A variety of trim calculation operations are achieved when 
this control location Is nonzero. If this control option Is 1. , the Glauert 
variable Inflow Induced velocity components Vq, Vj|_^, and v^g (locations 80-82) 
would be varied In addition to those parameters shown In Table II to effect 
momentum balances In thrust, pitching moment, and rolling moment as well as 
the specified trim. If the control option Is 2., only the velocity components 
would be varied (Instead of the control euigles, e A^_g, euid ^g)j but In 

an analogous manner. If the control option Is the trim calculation 
would again only vary the velocity components, but to achieve moment\im balance 
only, with no specified trim. 



APPENDIX III 


PROGRAM g400 OUTPUT DESCRIPTION 


The complete printed output generated hy the g 400 program can be 
classified into the following five ma^or categories; 

I. Listing of Input Data 

II. Parameters Calculated from the Input Data 

III. Results of Solution Part I - Eigensolutions 

IV. Results of Solution Part II - Time-History Solution 

V. Results of Solution Part III - Transient Spectral Stability Analysis 

This appendix will describe the pertinent output pages associated with each of 
these categories. It should be noted that while output will always be generated 
for the first two categories, output for the remaining categories depends upon 
the optional utilization of solution parts I, II, and/or III. The subsections 
which follow describe, in turn, the details of each of these five categories. 


Listing of Input Data 

Output in this category includes a listing of the static airfoil data (if 
any is inputted), a descriptive listing of the "inertia, elastic, geometric 
and other operational" (Loader Format) data, and a listing of the harmonic 
variable inflow: blocks I, II, and IV, respectively, of the input data des- 
cribed in Appendix II. If static airfoil data is Inputted, then a listing of 
this data will be outputted for c^, and <^mn/h with the format shown 

in Sample Page 1. where each column represents data at one Mach number. With- 
in each column the first line gives the number of angle-of-attack/ aerodynamic 
coefficient pairs defining the functionality; the second line is the Mach 
number, and the ensuing line pairs are the angle-of-attack/ aerodynamic coeffi- 
cient pairs, where the angles -of-attack are in degrees. This output closely 
follows the input format described in Appendix II. 

A description of the Loader Format data output is omitted herein since 
this output merely duplicates the description already given in even greater 
detail in Appendix II. In Sample Page 2 is shown the listings of the inputted 
harmonic variable inflow, both by haimonic components and by azimuthal varia- 
tion for each of the (maximum of) fifteen blade segments. As with the input 
format, the inflow is defined positive up, has the units of ft/sec and a 
conventional positive Fourier series representation is used. 
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Parameters Calculated from the Input Data 


Sample Pages 3 and If list a variety of intermediate parameter calcTilations 
which, for the most part are directly applicable only to the composite " 

bearingless rotor. Sample Page 3 presents the matrices and vectors used to 
define, respectively, the elastic bending characteristics of the torq.ue tube 
and the inboard (snubber) end vertical deflection of the torque tiibe. For 
each, a quadratic variation with pitch angle is assumed so that the total 
bending stiffness matrix, FDEFL, and inboard deflection vector, Z5I» are 
formed in the indicated manner wherein the total, blade pitch angle at the 
75 percent span, TTHO, is taken in radians. Furthermore, all the results out- 
putted for these two quantities are nondimens ional (see the superscripts sec- 
tion of the List of Symbols for the proper nondimens ionalization) . The matrix 
FDEFL duplicates matrix S, whose rows and columns are defined by eqmtion 
( 108 ) in the text; vector Z 5 I shares the same columnar dependence with matrix 
FDEFL. Finally, at the bottom of Sample Page 3 are evaluations of FDEFL and 
Z5I, xjsing the quadratic representations, at a pitch angle equal to the 
inputted collective angle, 9 . 75 R. It should be noted that the functionaliza- 
tions of these quantities to quadratic form is accomplished by matching 
exactly the quantities calctilated with pitch angles equal to the collective 
angles. 

The first group of output parameters given in Sample Page 4 consists of 
the (nondimensionaJ.) coefficients defining the nonlinear torsional stiffness 
characteristics of the flexbeam; they correspond directly to the coefficients 
given in equation (l30a). Again the standard nondimensionalization is used. 

The second group of output parameters on the sample page consists of the 
quadratic functionality coefficients for the pitch-flat and pitch-edge 
couplings and the angular motions of the torque tube relative to the flex- 
beam at the juncture for each of the selected flatwise and edgewise modes. 

The quadratic functionality on pitch angle duplicates that used for FDEFL and 
Z5I described in the above paragraph. Note that RELMN and RELMV are calcu- 
lated only on the basis of an infinitely rigid torque tube and, hence, are 
omitted when the redundant analysis option is invoked. The third group of 
output parameters consists of the specialized elastic description of the 
inboard end of a cantilevered torque tube (see Figure 6 and equation ( 83 ) )• 

The outputted spanwise distributions are the static flatwise deflection shape 
(pseudo- flatwise mode) and its two (nondimehs ional) spanwise derivatives. 

The fourth group of output parameters are the effective torsional springs 
K©i, K 02 > and defined by equations^( 63 b), ( 88 ), and ( 85 ), respectively, 

and the flexbeam stiffness , GJ and Tk^ , respectively. The units of the 
springs are Ib-ft/rad and those of the stiffnesses are Ib-ft^. 
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FLCX-BEAe TORSIONAL STIFFNESS CHARACTERISTICS 
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In Sample Page 5 ai*e shown typical modal information for the inputted 
flatwise and edgewise bending modes. For each such mode the ( nondimens ional) 
modal frequency, pitch- flat (or pitch-edge) coupling and inputted mode shape 
and spanwise derivatives are listed. In addition the listing presents the 
derived incremental deflection vectors which accoimt for blade twist (see 
eq,uations (7) thro\ij^ (10)), Within a flatwise modal information group, the 
DVB and DVE arrays correspond to those first order spanwise functions due 
to built-in twist and torsional modal twist, respectively. Ihe DWWBB, DWWBC 
and EWWCC arrays correspond to those second order ^W functions due to the 
combinations of built-in twist with itself, built-in twist with pseudo-torsion 
mode (control) twist, and pseudo- torsion mode twist with itself, respectively. 
The various arrays, DV2BP, DV2EP, DWW2BBP, DWW2BCP, and DWW2CCP are the first 
spanwise derivatives of the second components (those with superscript "2") of 
the above discTOSsed arrays, DVB, DVE, DWWBB, DWWBC, and DWWCC, respectively. 
Die UWE nonlinear deflection arrays correspond to the bracketed integral 
function defined in equation (36). Within the edgewise modal information 
group, the various arrays, DWB, DWE, DVBB, DWBC, and DWCC etc., correspond 
to similarly defined spanwise functions involving twist and the edgewise 
modal deflection and spanwise derivative arrays. 

The first group of output parameters presented in Sample Page 6 consists 
of the inputted torsion modal arrays together with the derived pseudo-torsion 
mode shape (as defined by Figure it) and spanwise derivative. The remaining 
output parameter group consists of the spanwise distributions of various 
pertinent aeroelastic quantities. The X and XDEN arrays are the nondimen- 
sional distances of the centers of the segments from the offset and rotor 
axis, respectively. The units of the CHORD array are feet. The angle of 
attack descriptors TEIETA-AERO, PEH and ALPflA are, respectively the geometric 
aerodynamic pitch angle, the inflow angle, and the resulting section angle- 
of-attack, all in degrees. These angles are calculated using the inputted 
initial conditions on azimuth angle and on the response variable deflections 
and velocities, Bie resulting Mach number and aerodynamic coefficients are 
used to define the perturbational airloads used in the eigensolution. The 
qiantlty KAPPA/u is the spanwise variation in aerodynamic moment damping 
coefficient which when multiplied by the local pitch rate approximates the 
potential flow unsteady pitching moment coefficient. The quantity (y10c/4)/c 
is the spanwise distribution of quarter chord offset from the reference axis 
nondlmenslonallzed by chord . The dynamic and structural quantities in the 
third group of this output page includes the QUAD array which constitutes the 
integration weighting numbers for spanwise integration. The THETA-STR array 
is the pitch angle distribution of the structural principle axes and has the 
units of degrees and, in general, differs from the aerodynamic pitch angle 
distribution. The two arrays, TWIST-BLT and TWIST-TOT, are the nondimensional 
structural twist rate distributions of the built-in twist and the total twist 
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(including elastic respcx^e and control inputs), respectively;, tbese arrays 
have the units of radians. The quantities TENSB, EIYB, Tgrzu and MASSB are, 
respective^, the blade tension, flatwise bending stlffhess, edgewise banding 
stlffhess and nass distributions all nondhnenslonallzed in the standard 
sense. Oiie (yiONA)/c and (yiOCG)/C arrays are, respectively, the edgairise' 
bending neutral axis and mass center offset distributions ^ 

by chord. 
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Results of Solution Part I - Eigensolutions 

Sainple Pages 7» and 9 present the pertinent details of the 
eigensolutions. Ihere are three distinct eigensolutions calculated etnd the 
three results are headed, respectively, by the following titles: 

1. 031UIICATED (LUJEAR TEEl® ONLY) VACUUM CASE 

2. LINEARIZED NONLINEAR VACUUM CASE 

3. LINEARIZED NONLINEAR NONVACUUM CASE 

In the first case, all terms nonlinear in the response variables and 
terms involving aerodynamic loadings are omitted. In the second case, those 
nonlinear inertial and elastic terms omitted in the first eigensolution are 
retained and linearized about the inputted response variables (initial condi- 
tions on response variable deflections and velocities — locations 67 throu^ 

70 and 280 throTo^ 303). Note that if these locations are inputted equal to 
zero this second eigensolution is automatically omitted. The last case 
includes the locally linearized perturbational airloads with the linearized 
formulation of the second eigensolution. As shown in Sample Pa^e 7j the 
beginning of each eigensolution lists the A, B, and C matrices, which pre- 
multiply the acceleration, velocity and displaceanent response vectors, res- 
pectively, to define the eigenproblem. The extracted eigensolutions (coupled 
roots and frequencies) are then listed by root pairs. At the bottom of 
Sample Page 7 and the top of Sample Page 8 eure presented the typical format 
for the case where each of the root pair is real. For each root the eigen- 
vector or GENERALIZED MODE SHAPE is presented, normalized to the Isurgest 
amplitude. The number of elements to this vector is identical with the dimen- 
sion of the A, B, and C matrices, and represent, in order, the coupled rela- 
tive responses of the flatwise, edgewise, and torsion modes. The arrays labeled 
mysiCAL M)DE SHAPE consist of the relative spanwise distributions of inplane 
(Y5), and out-pf -plane (Z5) and pitching (THETA) components of the coupled 
mode shape. The Y? and Z^ deflections are nondimens ionalized by blade radius, 
whereas THETA is in radians. 
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PART I. f iGiKSOLUTlONs OF VARIOUS LINEARIZATIONS OF EQUATION SET - CMAPACTER IS T1 C ROOTS AND COUPLED HOOC SHAPES 


truncated ilinear terms onlyi vacuum Case 


A matrix 
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Sample Page 7 
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At the hottom of Sample ‘ Peige 8 is presented the typical format for the 
case where the root pair consists of complex conjugates. For this case, 
in addition to the real and imaginary components of the roots, the eq.iiivalent 
critical damping ratio, ZETA, and undamped natural frequency, WN, are given. 
Note that for a complex pair of roots, the eigenvector is also complex; 
the sign on the imaginary component corresponds to the root with the 
positive imaginary part. Similarly, the physical mode shape is given. Which 
additionally contains the velocity component distribution to account 'for 
the generally nonuniform phasing along the span. The velocity consonants 
are nondimensionalized by tip speed. 
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Shoiild one of the roots in the nonvacuum eigensolutlon he unstable, as 
Indicated by a positive root or real part of a complejc pair, an output listing 
of the force phasing matrices appropriate to the instability is generated and 
outputted Eis depicted in Saii^>le Page 9* ^ese matrices, having the same 
size as the A, B, and C dynamic matrices, enable the various destabilizing 
forces to be identified; descriptive material for their definition and inter- 
pretation are contained in Reference 13. 


A PHASING matrix 
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Results of Solution Part II — Time History Solutions 

Sample Page 10 presents a variety of response and load quantities 
defining the transient aeroelastic responses. The first row of parameters 
following the page title presents, for the subsequent time-histoiy solution, 
a listing of the parameters defining the flight condition. These parameters 
consist of the various control angles (in degrees), the inflow and advance 
ratios, and the nondimensionaJ.ized values of the "momentum" induced velocity 
components. The remainder of Sample Page 10 comprises the typicail azimuthal 
listing; this listing is outputted for every azimuth angle which is a multiple 
of the print azimuth increment, input item no. 23. 
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The first of the four groups of result quantities on this sample page 
lists the spanwise distributions of the pertinent aerodynamic quantities. 

The inflow euid total section angles of attack PHI and ALPHA, respectively., are 
in degrees. The MACH NO., CL, CD, and CM are self-explanatory and nondimen- 
sional. The airload distributions in the Z5 and Y5 directions, SAZ5 and SAY5, 
respectively, have the units of Ib/in., the aerodynamic pitching moment 
distribution, MAX5, has the uriits of in.-lb/in. The quantities SDZ5, 

SDY5, and MDX5 are "semi-dynamic" load distributions. These distributions 
are dimensionally similar to those above described resulting from aerodynamics, 
but instead arise from all the dynamic effects except the doubly time differentia - 
ated ones (see equations Ulc, 4lb and li4a, respectively). The quant it^''MEX9 
is the elastic torsion moment distribution which consists of those torsion 
couplings arising from AEI, the tension-neutral axis offset and other twist 
related elastic effects; it too, has the iihits of in.-lb/in. 

The second group of result quantities consists of the instantaneous 
generalized excitations, XI, and the generalized accelerations, velocities, 
euid displacements (Q**,Q*, and Q, respectively) for each of the modal response 
variables selected, all appropriately nondimensionalized. The generalized 
excitations are defined to be the elements of the right hand side of equation 
(3U). 

’The third group of result quantities are comprosed mainly of blade 
deflection and stress distributions. The vertical and inplane deflections 
are those in the Z5 and Y5 directions, respectively. All stress quantities 
have the units of lb/in.2, whereas the torsion moment has the units of lb-in. 

It should be noted that, over the flexbeam-torque tube span, the flatwise 
and edgewise stresses outputted are those only for the flexbeam whereas the 
torsion moments and stresses outputted are those only for the torque tube. 

The last two columns are the distributions of the A and B parameters needed 
to define the unsteady airloads (see Reference 9 ). The fourth group of 
result quantities at the bottom of the sample page consists of miscellaneous 
deflection, load and stress results for the flexbeam and push-rod. 

The line titled SPAR/FLEXURE PARAMETERS presents similar stresses and 
torsion moments for the flexbeam at the indicated spanwise location; note, 
however, that the outputted torsion moment is for the flexbeam immediately 
inboard of the Juncture. Of the remaining four quantities, PUSH-ROD (RELATIVE) 
DEFL. (in.), TORQUE_TUBE ROOT DEFL. (in.), and TORQUE TUBE ROOT SHEAR (lb.) 
pertain only to cantilevered torque tube configurations wherein the "wobble 
mode" option is invoked (input location 985). All three quantities are defined 
in the positive Z5 direction. The quantity PUSH-ROD LOAD (lb) is the upward 
(+Z5) directed load the push-rod exerts on the pitch horn/push-rod attachment 
point; this qusuatity is csJ.culated for all blade types. For composite bear- 
infeless rotors, the push-rod load accounts for the total blade torsion moment 
at the root less that torsion moment resisted by the flexbeam. 
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PART IT. HISTORY SOLUTION OF CONPLLTL INONLINFAR) LCUATION SCT - Ae*OCLASTJC TRAKSKMT RCSPONSLS 
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After the time-history solution has either converged to periodicity or 
run to maximum flapping tfials (input location 10) vario'us integrated loads 
are calculated for one final hlade revolution to form the aerodynamic perfor- 
mance and stress results depicted in Sample Page 11. For each of eight (8) 
performance quantities results are presented in nondlmenslonal coefficient 
form, in nondlmenslonal form divided by solidity, and in actual dimensional, 
form. Note that ten (10) dimensional quantities are listed and the units are 
lb for forces and Ib-ft for moments, as appropriate. Uie quantity EQU. DRAG 
(lb) represents the combined power expended by the rotor due to rotor rotation 
(torque) and translation (drag) divided by flight speed. 

GSie next line duplicates the parameters defining the flight condition and 
includes four (4) additional quantities which depend on the integrated perfor- 
mance for evaluation. At the beginning of the time-history calculation it is 
not known which part of the inflow ratio being •used is due to ram effects and 
which due to momentum Induced effects. Once the integrated rotor thrust is 
calculated, however, the induced portion of the inflow can then be calculated 
using the simple usual momentum formula derived for flight in an infinite con- 
tinu'um (Reference 14). Ihe complementary portion of the Inflow represents the 
ram effect from which the shaft angle- of-attack ALPHA S, in degrees, can be 
calciilated. Ihe quantity VEL ACT. is the actual forward flight velocity, in 
knots, consistent ■with the ad'vance ratio used and the shaft angle of attack. 
For finite forward flight speeds EQU. l/D is the lift di-vlded by the equi"va- 
lent drag; for hovering cases this q'uantity is the figure of merit. PAR. 

AREA, the rotor parasite (drag) area, in square feet, is the rotor drag di-vlded 
by dynamic pressure. The line titled CORRECTIQETS DUE TO WIND TUNNEL WALL 
INTERFERENCE; consists of recalculations of those qxxantities which depend on 
the Induced portion of the inflow wherein the induced inflow is calculated 
using the formulae derived for flight in wind tunnels or in ground effect 
(Reference 11). The remainder of Sample Page 11 consists of reductions of the 
various stresses (Ib/in.^) and the push-rod load (lb.) to median and •§■ peak- 
to-peak values over the final rotor re-volution. 
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AEPQOYNAHK PlRFORNANCE AND STRESSES 
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Should ma;)or (trim) Iterations he used (see description of Input Items 
6o and 62 through 6? In Appendix II) output depicted on Sample Page 12 irlll 
he generated hy the program. Ihe first line consists of the , zeroth, first 
cosine and first and second sine harmonics of first flatwise mode response, in 
radians, and an estimate of an effective angle-of>attach on the retreating 
hlade side (i)r > 270°), In degrees. Ihe nonzero elements of the depicted (G) 
MATRIX give, for each row, the partial derivatives of the four trim quantities 
9 rm» respectively) with respect to the four control quantities 
being used (0 «=» \g» and (sin a^), or v^^, Vj^^, v^^g and (sin ofg)), for 

each respective column. !0^e elements of this matrix are formulated using 
approximate linear strip theory with reversed flow effects and are calculated 
for either set of control quantities, as appropriate. Ohe ERROR VECTOR con«> 
slsts of the differences between the four requested tirlm quantities and those 
achieved in the preceedlng tlme>history. Ihe two lines depicted give the error 
vector In dimensional (lb and Ib-ft) and nond Imenslonal fozms, respectively. 
!Ihe CORRECTION VECTOR consists of those changes to the control quantities 
which should null the above described error vector. iQie correction vector Is 
obtained from the premultlpllcatlon of the inverse of the G matrix with the 
error vector, but the corrections are scaled, if necessary, to prevent control 
changes of more than 2 degrees within any one Iteration. Ohe control param- 
eters whose increments are depicted in this output page are, in respective 
order: 9 Aj^g, Bj^g, (sin ofg), X, C^, v^, Vj^^, and v^^; the first three have 

units of degrees and the remainder are dimensionless or nondimenslonallzed. 


IkO 



trim XTERATIONi 
U) MATRIX 

> OHIO s 

X.XXXE-YY 

ONlC = 

XoXXXE-YY 

owls : X.XXXE-YY 

MIS s X.XXXC-TT «(t70r a XX.XX 

X.XXXE-tt 

X.XXXE'YV 

X.XXXE-YY 

x.xxxe-yy 

0.000 

0.000 0.000 

0.0«0 

0.000 

x.xxxe-tt 

X.XXXE'VY 

X.XXXE-YY 

x.xxxe-yy 

0.000 

0.000 0.000 

0.000 

0.000 

X.XXXE-TT 

X.XXXE-YY 

x.xxxe-yy 

x.xxxe-yy 

ObOOO 

0.000 0.000 

0.000 

0.000 

X.XXXE-TV 

XfXXXE-YY 

X.XXXE-YY 

x.xxxe-yy 

0.000 

0.000 0.000 

0.000 

0.000 

ERROR VECTOR 








X.XXXE-TY 

X.XXXE-YY 

X.XXXE-YY 

x.xxxe-yy 

0.000 

0.000 0.000 

0.000 

0.000 

x*xxxe*YY 

X.XXXE-YY 

x.xxxe-yy 

x.xxxe-yy 

0.000 

0.000 0.000 

0.000 

0.000 

COMKCTtON VECTM 







x.xxxe*w 

XeXXXE*YT 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

X.XXXC-VV X.XXXC-VT X.XXXE-TT 

XeXXXC-YY 
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Once the time-history solution has converged to periodicity and aJi major 
iterations have been completed, the program optionally performs harmonic anal- 
yses of the azimuthal variations of various response quantities (see descrip- 
tions of input items 48 and 73 through 75 in Appendix II ). Ihe outputs of 
these harmonic analyses are depicted in Sample Pages 13 throxigh 15. In each 
of these sample pages the harmonic information for each response variable is 
contained in the appropriate horizontal band of five rows, llie harmonics are 
listed by columns up to a maximum of 10 harmonics. All harmonic analysis out- 
put depicted on these sample pages assume a negative harmonic content form in 
keeping with the (negative) harmonic form conventionally as^-amed for the blade 
pitch control and rigid flapping angles. For each harmonic of response vari- 
able five quantities are outputted; these quantities are, respectively, the 
cosine and sine components, the equivalent amplitude and phase angle, and 
lastly, the amplitude of the harmonic relative to all the other harmonic ampli- 
tudes outputted. Sample Page 13 depicts the harmonic analyses of the dimen- 
sionless modal response variables selected wherein QW(l), QV(k) and QT(j) are, 
respectively, the (l) flatwise, (k) edgewise and (j) torsional uncoupled mode 
responses. 

Sample Page l4 depicts the harmonic analyses of the total shears and 
moments exerted by one blade to the hub. In contrast to the steady hub loads 
listed in the AERODYNAMIC PEEPOEMANGE AND STRESSES output (Sample Page 11 ) 
which are calculated by integrating only the aerodynamic load distributions, 
the total hub loads which are herein harmonically analysed are calculated by 
similarly integrating the combined aerodynamic and the dynamic load distribu- 
tions. Bie longitudinal, lateral and vertical hub shears comprising the first 
three quantities of this sample page all have the dimensions of lb and are 
defined in the X 2 “(aft), y^- (starboard), and z^-(up and along axis of rotation) 
axis directions, respectively. The roll, pitch and yaw moments comprising the 
latter three quantities on this sample page have the dimensions of Ib-ft and 
are defined positive (using the right-hand rule) about the x.-, y^^-, and z^- 
axes, respectively. Note that the aerodynamic rolling moment whose output is 
depicted in Sample Page 11 is defined positive starboard side down and is 
opposite from the harmonically analyzed total rolling moment depicted in Sam- 
ple Page l4. Sample Page 15 depicts the harmonic analysis of the flatwise 
stresses at the center of each of the spanwise segments. A similar output 
listing is provided for both edgewise and torsional stresses. 
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HARHONIC ANALYSIS OF BLADE RESPONSES 


QUl 


QW2 


OU3 


QVl 


QV2 


QTl 


AO 

1 

2 

3 

4 

5 

b 

7 / 

8 

9 

10 

A.X5«XE-YY 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX .X 
.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX .X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 

x.xxxE-vr 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX. 

X.XXXE-YY 
X.XXXE-YY 
X.XXXF-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X ' 

• XXX 

x.xxxe-yy 

X.XXXE-YY 

x.xxxE-yy 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
X XX .X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
. .XXX 

X.XXXE-YY 

X.XXXE-YY 

x.xxxE-yy 

XXX.X 
• XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxe-yy 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 

X.XXXE-YY 

x.xxxe-yy 

XXX.X 
• XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX .X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxe-yy 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

x.xxxe-yy 

X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY. 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxE-yy 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.’XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X .XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X .XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxe-yy 

X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxe-yy 

X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X ,xxxe-yy 

X.XXXE-YY 

X.XXXE-YY 

X.XXXE-YY 

XXX.X 

.XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X .X XXL -YY 
X.XXXL-YY 
X .XXXL-YY 
XXX.X 
• XXX 

X .XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X .XXXE-YV 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

x.xxxe-yy 

X, XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X ,XX XE-YY 
X .XX XE-YY 
X.XXXE-YY 
XXX.X 
. XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-YY 
X.XXXE-YY 
X.XXXE-YY 
XXX.X 
• XXX 

X.XXXE-Y 
X.XXXE-Y 
X.XXXE-Y 
XXX.X 
• XXX 
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H*»MOkIt ANALYSIS OF HU« SHlANS AND nOMCnTS 



AO 

1 

2 

3 

tl 

s 

* 

7 

II 

♦ 

10 

L0N6. SHR 

xx.x 

xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 



xxxxx;x 

xxx*x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 



xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 



xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 



x.xx 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

.XX 

• XX 


LAT. SHR 

xx.x xxxxx.x 

xxx.x 

xxx.x 

XX.X 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 


xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 


xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 


xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 


x.xx 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 


VCRT. 

SHR 

xx.x 

xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 




x.xx 

• XX 

• XX 

.XX 

• XX 

• XX 

• XX 

• XX 

• XX 

.XX 

ROLL 

HOHT 

xx.x 

xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

XX. X 

xx.x 

xx.x 




xxxxx.x 

xxx.x 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

rxx.xx 

xxx.xx 

xxx.xx 




x.xx 

• XX 

• XX 

• XX 

• XX 

.XX 

• XX 

• XX 

• XX 

• XX 

pitch 

mohT 

xx.x 

X xxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

*xx.x 

xx.x 

xx.x 

jlX.X 

AX.X 




xxxxx.x 

xxx.x 

xxx.x 

xx.x 

XX.X 

xx.x 

xx.x 

xx.x 

xx.x 

AV.X 




xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 




x.xx 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX 

.XX 

• XX 

• XX 


7AW 

HOHT 

xx.x xxxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

xx.x. 

xx.x 

xx.x 

xx.x 



X xxxx.x 

xxx.x 

xxx.-x 

xx.x ‘ 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 



X xxxx.x 

xxx.x 

xxx.x 

xx.x 

xx.x 

xx.x 

AX .y . 

xx.x 

xx.x 

xx.x 



xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

xxx.xx 

■ xxx.xx 



x.xx 

• XX 

• XX 

• XX 

• XX 

• XX 

• XX . 

• XX 

• XX 

• XX 
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HIR»*ONIC ANALYSIS OF FLATtflSC STRESSES 


N 

X CEN 

AO 

1 

? 

3 

k 

s 

6 

7 

8 

9 

10 

1 

• xxxx 

xxxx.x 

xxxx.x 

XXXX.X 

xxxx.x 

xx.x 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

xxxx.x 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




XX. X 

xx.x 

xx.x 

xx.x 

xx.x 

XX.X 

xx.x 

xx.x 

xx.x 

xx.x 




• XXX 

x.xxx 

.xxx 

• xxx 

• xxx 

• xxx 

• xxx 

•xxx 

• xxx 

•xxx 

2 

• xxxx 

xxxx.x 

xxxx.x 

xxxx.x 

xxxx.x 

xx.x 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

xxxx.x 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

XX.X 

xx.x 

xx.x 

xx.x 

xx.x 




.XXX 

x.xxx 

• xxx 

• xxx 

• xxx 

• XXX 

• xxx 

• xxx 

• xxx 

• xxx 

3 

.XXXX 

xxxx.x 

xxxx.x 

xxxx.x 

xxxx.x 

xx.x 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

xxxx.x 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




XX. X 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




.XXX 

x.xxx 

• xxx 

.xxx 

• xxx 

• xxx 

• xxx 

.xxx 

• xxx 

.xxx 

« 

• XXXX 

xxxx^x - 

xxxx.x 

xxxx.x 

xxxx.x 

xx.x 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

xxxx.x 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

XK.X 




xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




• XXX 

X.XXX 

• xxx 

• xxx 

• xxx 

• xxx 

• xxx 

.xxx 

• xxx 

• xxx 

5 

• xxxx 

xxxx.x 

xxxx.x 

xxxx.x 

xxxx .X 

xx.x 

xxx^x 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

xxxx.x 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

XX.X 

xx.x 

xx.x 

xx.x 

xx.x 




• XXX 

x.xxx 

• xxx 

• XXX 

• xxx 

• xxx 

• xxx 

• xxx 

• xxx 

• xxx 

6 

• XXXX 

xxxx.x 

xxxx.x 

xxxx .X 

xxxx.x 

xx.x 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xxxx.x 

XXXX.X 

xxxx.x 

XXX.X 

XXX.X 

XXX.X 

xx.x 

xx.x 

xx.x 

xx.x 




xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 

xx.x 




• XXX 

x.xxx 

• xxx 

.xxx 

• xxx 

xxx 

• xxx 

• xxx 

• xxx 

•xxx 


H 
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Results of Solution Part ill - 
Transient Spectral Stability Aaalysis 

Transient time-history solutions are often difficult to interpret for 
quantitative stability information. Ihis is due to the fact that the total 
responses so calculated inherently consist of several con^nent modes simulta- 
neously and transiently approaching (or departing from) multi-harmonic perio- 
dicity with a wide range of natural frequencies and inherent damping levels, 
!Ihe extraction of the component responses at discrete frequencies in order to 
examine their individual attenuation characteristics is the purpose of the 
transient Spectral Stability Analysis (tSSA) portion of Erogram g 400, Ihe 
details of this analysis, which utilizes Fourier Transform techniques, are 
beyond the scope of this report but are treated in References 15. 

Essentially the TSSA first performs Fourier transformations of selected 
time-history data strings, which have been previously generated in the time- 
histoiy solution portion of the analysis (Solution Bart II) and appropriately 
saved. The purpose of the Fourier Transform is to identify, within these 
time-histories, those frequencies whose amplitudes are relatively largest and 
and which are herein denoted as "resonances". The TSSA then calculates the 
transient behavior of the extracted amplitudes of these resonances over the 
time-history time interval and estimates equivalent linear stability indices 
(characteristic exponent, critical damping ratio, and time to half amplitude). 

Sample Rages l6 through l8 depict the output typically generated by the 
TSSA. The sequence of output depicted is duplicated for each of the transient 
response channels selected (see input locations 84 through 86 , Appendix II). 
Sample Page l6 depicts the output generated by the Fourier Transform frequency 
identification portion of the TSSA. Shown at the top of the page is the tran- 
sient response channel being analyzed and the frequency range wherein resonance 
identification is desired (input locations 8? and 88). The series of five out- 
put items to follow consist of parameters defining the nimierical Foiirier Trans- 
form; note that the results of the TSSA incoiporate a time nondlmensionaliza- 
tion based on rotor speed, fi. The tabulation of the Fourier Transform follows 
wherein, for each frequency (harmonic of the fundamental as determined by the 
total nondimensional time interval), the real and imaginary parts, the square 
of the azzplitude and the logarithm to the base 10 of the amplitude eire out- 
putted. Generally, this tabulation will consime more them the one page indi- 
cated in Sample Page l6. After this listing is completed, those frequencies 
and their respective square amplitudes which are foxuid to be resonances, as 
defined above, are listed. 
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PIUT 111. TRAMSICNT SPECTRAL STABILITY ANALYSES OF SELECTED AEROEL^STIC TRANSIENT RESPONSES 


TRANSIENT RESPONSE CHANNEL NO. XX - FOURIER TRANSFORM AND RESONANT FREOUENCV IDENTIFICATION 

OFSIREO FREQUENCY RANGE = .XX TO X.XX < /REVi 


NUMBER OF POINTS IN TIME SERIES 
INUI TIME INCREMENT 
LENGTH OF INTERVAL 

fundamental frequency 
highest frequency 


XXX 

•XXXXX RAO 
XX.XXXXX RAO 
.XXXXX /REV 
XX.XXXXX /REV 


TABULATION OF FOURIER TRANSFORM BY HARMONIC OF FUNDAMENTAL FREQUENCY 







TRANSFONM 

LOG TO BASE 1C 

HAR*»0NIC FRE&UFNCY 

COMPLFx FOURIER TRANSFORM 

magnitude 

OF magnitude 


Q 

.xxxxxx 

X. XXXXXX XE-YY 

•OODOOOD I 

X.XXXXXXXE-YY 

X .XXXXXXXE-YY 


1 

•xxxxxx 

X. XXXXXXXF-YV 

X.XXXXXXXt-YYl 

X.XXXXXXXE-YV 
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Aside from their several nonlinearities, the dynamic equations of motion 
of helicopter rotor blades implicitly contain several linear terms with period 
coefficients, which arise from the periodic character of the airloads in for- 
ward flight. It is not unexpected then, that the aeroelastic time-history 
responses generated by these equations should manifest Ploquet Theory charac- 
teristics (see Reference l6). In particular, the Fourier Transform is capable 
of identifying '*multiple resonances” which are separated hy (plvis or minus) 
multiples of the rotor frequency and which would be found to have approximately 
the same damping level. Hence, the resonant frequencies fotuid by the resonance 
search are further screened to extract only those frequencies with distinct 
nonlnteger values and which, within the set having the same noninteger values, 
have the largest transform magnitudes. These exd;racted frequencies are herein 
denoted ’’fundamental resonances” and are the only ones examined further for 
stability in the TSSA. 
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Saaiple page 17 depicts the resvilts of frequency fine tuning and response 
stability estimation for each of the fundamental resonances extracted earlier 
in the ISSA. Ihe results for each of these freq,uencies are presented in 
columnar fashion. Ihe top horizontal blocks of output represent the fre- 
quency fine-tuning results. Of most practical importance are the values 
labeled OPTIMIZED PREQUENCT which are, in nondimensional (per rotor rev) form, 
the best estimates of the frequency of the fundamental resonant frequencies. 
Dxese frequencies are obtained by an optimization technique, the details of 
tdilch are beyond the scope of this report. The remainder of the output 
depicted on this sample page (for each fundamental resonance) consists of 
three horizontal blocks of output representing various estimates of the effec- 
tive damping characteristics. Ihese three types of blocks are best explained 
by first describing Sample Page l8^ TMs sample page depicts, columnarly for 
each of the fundamental resonances indicated in Sample Page 17, the natural 
logarithm of the magnitude of resonant frequency content at each (nondimen- 
sional) time indicated. If these amplitude logarithms attenuate with time, 
then that frequency content (mode) is deemed stable, and conversely the slope 
of that attenuation with time is a measure of the effective linear damping; in 
the analysis this slope is obtained by a simple least-square fit. It may 
happen that the variation of amplitude logarithms with time is neither mono- 
tonic increasing or decreasing in which case a condition of maximum or minitmm 
amplitude is defined. Ey weighting the least-square fit either uniformly or 
with an appropriate function accentuating the initial or terminal ends of the 
amplitude logarithms data string, the three latter horizontal blocks of output 
depicted in Sample Page 17 are generated. Within each of these blocks, the 
first quantity depicted is the nondimensional CHAEACTERISTIC EXPONENT, which 
is analogous to and interpreted in the same way as the real part of the eigen- 
value discussed in the output for Solution Part I. The REVS to (MAX/MIN) AMPL, 
is an indication of the asymptotic behavior of the component response. STAN- 
DARD DEVIATION is the root-mean- squared error achieved in the least-square 
curve-fit and is an indication of the regularity of the amplitude logarithm 
function depicted in Sample Page l8, and of the accuracy of the stability 
estimation. Based tpon the OPTIMIZED FREQUENCY outputted at the top of the 
sample page, the equivalent CRITICAL DAMPING RATIO is calculated from the 
characteristic exponent using standard formulae. Finally, the output item 
labeled REVS TO HAIF AMPLITUDE is the third alternate way in which the equiva- 
lent linear damping result is presented. 

m Sample Page 19 is depicted the typical additional page of output 
generated at the beginning of every case following the first case of a multiple 
case sum. The two coltimns depict, respectively, the location numbers and data 
values for the newly inputted data distinguishing the present case from the 
previous one. This feature is intended solely as an ease of usage output to 
assist in data management. 
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